Unit 11 


Eigenvalues and eigenvectors 


Introduction 


Introduction 


This unit develops the study of geometric transformations that you began 
in Unit 6. Much of our attention will be on finding the position vectors of 
points that are mapped by a linear transformation to scalar multiples of 
themselves. Here a scalar multiple of a point (x,y) is another point of the 
form (kx, ky), for some real number k. 


Consider, for example, the linear transformation 
r()=(4 8) @): 
Observe that 
(1 3) G)=@)-2(). 
1 3 1 2 1 
This computation shows that f maps the point (1,1) with position 


1 : è ; : i ; 
vector C) to a scalar multiple of itself. More precisely, it scales this point 


and the corresponding position vector by the factor 2, as illustrated in 
Figure 1. 


(2,2) 


Figure 1 A transformation mapping (1,1) to (2,2) 


A non-zero vector that is scaled by a linear transformation in this sort of 
way is called an eigenvector of the matrix representing the tranformation, 
and the scale factor is called an eigenvalue of the matrix corresponding to 
this eigenvector. Here you’ll learn how to find the eigenvalues and 
eigenvectors of 2 x 2 matrices. In fact, you’ll see that finding the 
eigenvalues of a 2 x 2 matrix comes down to solving a quadratic equation, 
and finding eigenvectors corresponding to these eigenvalues is a matter of 
solving pairs of simultaneous equations. 


Later in the unit you'll get a taste of some of the uses of eigenvalues and 
eigenvectors in modelling real-world situations. For example, you’ll see 
how the eigenvalues and eigenvectors of matrices can be used in an efficient 
method for calculating powers of matrices, which has applications to 


97 


Unit 11 Eigenvalues and eigenvectors 


98 


modelling the sizes of human or wildlife populations in a subject called 
population dynamics. You'll also find out how the eigenvalues and 
eigenvectors of a matrix are employed in solving systems of differential 
equations; systems that have important applications to the biological and 
physical sciences. 


1 Eigenvalues and eigenvectors of 
matrices 


In this section you'll find out what eigenvalues and eigenvectors of 
matrices are, and learn strategies for finding them. 


1.1 What are eigenvalues and eigenvectors? 


In Unit 6 you learned that a linear transformation of the plane R? is a 
transformation f whose rule is of the form 


f(x,y) = (ax + by, cx + dy), 


where a, b, c and d are real constants. You also learned how to represent a 
linear transformation by a 2 x 2 matrix, as follows. Let 


a(i i) a (3) 


The vector x is the position vector of the point (x,y), in column vector 
form. Since 


_ fa b\ (z\ _ fax+by 
axs (e a) (i) (era): 


it follows that the rule for f can be expressed in the column vector form 


f(x) = Ax. 


As you saw in Unit 6, one-to-one linear transformations map lines to lines. 
You didn’t learn about the images of lines under linear transformations 
that are not one-to-one in Unit 6, but in fact such transformations map 
lines to either lines or single points, though we won’t prove that here. We 
say that a line is an invariant line of a linear transformation if its image 
under the transformation is the original line, or just a point on the original 
line. In this unit we’re particularly interested in finding invariant lines of 
linear transformations that pass through the origin. Since linear 
transformations fix the origin, a linear transformation that has an 
invariant line that passes through the origin maps that line either to itself, 
or to the origin. 


For example, consider the linear transformation f represented by the 
matrix 


"i 
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You'll recall from Unit 6 that this transformation is called a (2,3)-scaling 
because it scales horizontal distances by the factor 2 and scales vertical 
distances by the factor 3. In particular, it fixes the x-axis and the y-axis, 
so both these lines are invariant lines of the transformation, and both lines 
pass through the origin. 


Since the z-axis is an invariant line of the (2,3)-scaling f, any point on the 
x-axis is mapped to a scalar multiple of itself. Here are some examples, 
using matrix and column vector notation: 


(a) = (0) AO) = 6) aC) =): 


The blue arrows and points on the x-axes in Figure 2 illustrate the first of 
these equations. (The other arrows and points in Figure 2 will be discussed 
shortly.) 


Figure 2 The effect of the (2, 3)-scaling 


These examples indicate that f scales any point on the x-axis by the 
factor 2. To verify this, choose any point (k,0) on the x-axis, and observe 
that 


k 2k k 
a (a) = (8) =2(0), 
so f does indeed scale (k, 0) by the factor 2. 


The y-axis is also an invariant line of f, so points on the y-axis are also 
mapped to scalar multiples of themselves, as the following calculations 
demonstrate: 


ae A JE ra 


The purple arrows and points on the y-axes in Figure 2 illustrate the first 
of these equations. The examples indicate that f scales any point on the 
y-axis by the factor 3. You'll carry out a calculation to check this assertion 
in the next activity. 
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Figure 3 Reflection in the 
line y= x 


Figure 4 Images of points 
under the reflection in the 
line y= x 
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Activity 1 Investigating the (2, 3)-scaling 


Let A = i AF Show that if x = a for some number k, 
then Ax = 3x. 


To summarise, we’ve found that for the linear transformation f(x) = Ax, 

2 0 
where A = & 3 
position vector of a point on the z-axis, and the equation Ax = 3x is 
satisfied whenever x is the position vector of a point on the y-axis. The 
first equation tells us that the x-axis is an invariant line of the 
transformation f, and each point on the x-axis is scaled under the 
transformation by the factor 2. Likewise, the second equation tells us that 
the y-axis is an invariant line of f, and each point on the y-axis is scaled 
by the factor 3. 


i; the equation Ax = 2x is satisfied whenever x is the 


Let’s find a similar pair of equations and invariant lines for another linear 
transformation that you met in Unit 6, namely the linear transformation f 
represented by the matrix 


A=(j aE 


This transformation is the reflection in the line y = a, illustrated in 

Figure 3 by a triangle and its image. This line is an invariant line of the 
reflection; in fact, each point on the line is fixed by the transformation. For 
example, the points with position vectors 


(i) G) =» ) 


all lie on the line y = x, and 


1 1 3 3 —1 —1 
E A) ae 
The line y = —a, which is perpendicular to the line y = x, is also an 


invariant line of f that passes through the origin. The effect of f on points 
on these two lines is illustrated in Figure 4. 


Each point on the line y = —a is scaled by f by the factor —1. For 
example, the points with position vectors 


i) 2 


both lie on the line y = —a, and 


Oe eee Bore) ae aa) 
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In the next activity you are asked to perform calculations that confirm 
that f fixes each point on the line y = x, and f scales each point on the 
line y = —2 by the factor —1. 


Activity 2 /nvestigating the reflection in the line y = x 


0 1 
Let A = (; a 


(a) Show that if x = (; 


), for some number k, then Ax = x. 


(b) Show that if x = = 


J for some number k, then Ax = —x. 


In each of the examples you’ve seen so far there has been 
a 2 x 2 matrix A, a vector x and a number A such that 


Ax = Xx. 


More generally, we have the following definitions, which apply to square 
matrices of any size (not just 2 x 2 matrices). 


Eigenvalues and eigenvectors 
An eigenvalue of a square matrix A is a number such that 
Ax = Àx 


for some non-zero vector x. The vector x is called an eigenvector 
of A corresponding to À. 


For each eigenvalue ÀA of A, the equation Ax = Ax is called the 
eigenvector equation of A for the eigenvalue A. 


For example, if 


2 0 1 
i= ) and = 


then, as you have seen, Ax = 2x, so 2 is an eigenvalue of A, and x is an 
eigenvector of A corresponding to the eigenvalue 2. In fact, any non-zero 
scalar multiple of x is also an eigenvector of A corresponding to the 

eigenvalue 2. This is because if we scale x by the factor k then we obtain 


the vector GI and 


alo) =2(0): 


as you saw earlier. 
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More generally, suppose that x is an eigenvector corresponding to an 
eigenvalue À of some matrix A (so that Ax = Ax). Now choose any 
non-zero scalar multiple of x, say kx, for some non-zero real number k. 
Then 


A(kx) = k(Ax) = k(Ax) = A(kx), 


which shows that kx is also an eigenvector of A corresponding to A. These 
vectors kx are the position vectors of the points on a line through the 
origin. This line is an invariant line of the linear transformation 
represented by the matrix A. 


We record this important observation about scalar multiples of 
eigenvectors in the box below. 


Scalar multiples of eigenvectors 


If x is an eigenvector of a square matrix A corresponding to an 
eigenvalue À of A, then any non-zero scalar multiple of x is also an 
eigenvector of A corresponding to the eigenvalue A. 


You’ve just seen that there are infinitely many eigenvectors corresponding 
to a given eigenvalue. In contrast, given an eigenvector x of a matrix A, 
there can be only one number àA that satisfies the equation Ax = Ax. So 
for a given eigenvector there is only one corresponding eigenvalue. 


You'll notice from the definition of eigenvectors that an eigenvector is a 
non-zero vector. This means that an eigenvector cannot equal the zero 
vector 0, which in two dimensions is 


o= (9). 


The zero vector is excluded because the equation 
AO = \0, 


where 0 is the zero vector in n-dimensions, is true for any n x n matrix A 
and any number A, as both sides of the equation equal 0. 


Origin of the terms ‘eigenvalue’ and ‘eigenvector’ 


The word ‘eigen’ is German, meaning ‘own’ or ‘self’, and was brought 
into major use in the early twentieth century by the German 
mathematician David Hilbert (1862-1943). However, in English the 
term ‘proper value’ was used instead of ‘eigenvalue’ and remained in 
use until the 1960s. Indeed, in 1967 the Hungarian-born 
mathematician Paul Halmos finally admitted defeat, as you can see 
from the following extract from the preface to his text A Hilbert Space 
Problem Book. 


1 Eigenvalues and eigenvectors of matrices 


For many years I have battled for proper values, and against the 
one and a half times translated German-English hybrid that is 
often used to refer to them. I have now become convinced that 
the war is over, and eigenvalues have won it; in this book I use 
them. 


(Halmos, P. (1967) A Hilbert Space Problem Book, D. Van Nostrand 
Co., p. x) 


Example 1 Finding an eigenvalue corresponding to a given 
eigenvector 


—1\. : : 
Show that the vector ( J is an eigenvector of the matrix 


d g 
2 -1/)’ 
and find the corresponding eigenvalue. 


Solution 


@. You need to show that Ax = Ax for the given matrix A and 
vector x, and for some number A, which you must find. First 
calculate Ax. @ 


“E92 


@. Write the vector a in the form À E for some 


-2(-?). 


So Z is an eigenvector of the matrix corresponding to the 


number à. ® 


eigenvalue —2. 


Activity 3 Finding eigenvalues corresponding to given eigenvectors 


Show that the vectors (i) and (i) are eigenvectors of the matrix 


l 


and find the corresponding eigenvalues. 
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1.2 Finding eigenvalues 


In the previous subsection you saw that if 


2 0 
a=(0 3): 


then the eigenvector equation Ax = 2x is satisfied when x = 6 , and the 


0 
eigenvector equation Ax = 3x is satisfied when x = (2) . This shows that 
2 and 3 are eigenvalues of A. 


However, it’s often the case that you won’t know (or you won’t be able to 
guess) eigenvectors of a matrix in advance. Here you’ll learn a systematic 
way to find the eigenvalues of a 2 x 2 matrix directly, without first 
knowing what the eigenvectors are. You'll see that a 2 x 2 matrix has 
either one or two eigenvalues, which might be real numbers, or they might 
be complex numbers that are not real numbers. In the next subsection, 
you'll use the eigenvalues of a matrix to find corresponding eigenvectors. 
Finally, at the end of this section, we’ll discuss how to find eigenvalues and 
eigenvectors of square matrices of size larger than 2 x 2. 


The method we use to find the eigenvalues of a 2 x 2 matrix A is to obtain 
a quadratic equation whose solutions are the eigenvalues. To this end, let’s 
begin by writing the eigenvector equation 


Ax = Xx, 


for a given eigenvalue A in an alternative form. First subtract Ax from 
both sides to give 


Ax —Ax=0. (1) 
The right-hand side is the zero vector 0 = Ge 


It would now be helpful to factorise the left-hand side of this equation by 
taking out the common factor x. However, we can’t do this directly, as this 
would leave the factor ‘A — A’ which doesn’t make sense: we can’t subtract 
a number from a 2 x 2 matrix. 


Instead, we can make use of the 2 x 2 identity matrix 


(9) 


This matrix satisfies Ix = x, so we can write equation (1) as 
Ax — \Ix = 0. 


Since A and AI are both 2 x 2 matrices, we can use standard properties of 
matrices to write the left-hand side as (A — AI)x. This gives an 
alternative, equivalent form of the eigenvector equation. 
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Equivalent form of the eigenvector equation 
An equivalent form of the eigenvector equation Ax = Ax is 


(A — AI)x = 0. 


This equivalent form of the eigenvector equation is valid when A is 
an n x n matrix, for any positive integer n, provided that I is interpreted 
as the n x n identity matrix, and x and 0 are n-dimensional vectors. 


We can use the equivalent form of the eigenvector equation to explain how 
to find the eigenvalues of a 2 x 2 matrix. To demonstrate the process, let’s 
write 


a b 
at. 
and let’s suppose that is an eigenvalue of A with a corresponding 


eigenvector x = 


As A and AI are 2 x 2 matrices, the expression A — AI also represents 
a 2 x 2 matrix. For convenience, let’s denote this matrix by B, so that the 
equivalent form of the eigenvector equation for the eigenvalue A is 


Bx = 0. 
If B were an invertible matrix, then you could multiply both sides of this 
equation on the left by BT} to give 
B-'Bx = B~'0. 
Since B~!B =I, the resulting equation simplifies to give 
x=0. 
But this cannot be so, because x is an eigenvector of A, so it’s not the zero 


vector. It follows that B is not an invertible matrix after all, and so it has 
determinant zero. Now, 


eee k J - (ò : = Ce re 
so 
det B = (a — A)(d— A) — be = X? — (a + d)À + ad — be. 
Since det B = 0, we obtain the quadratic equation 
à? =(a+@A\+4 ad—be=0. 
This quadratic equation is called the characteristic equation of A. 


We’ve shown that any eigenvalue of A is a solution of the characteristic 
equation of A. Let’s now reverse this argument to show that any solution 
of the characteristic equation of A is an eigenvalue of A. Suppose then 
that u is a solution of the characteristic equation, so that 


det(A — pI) = p? — (a + d)u + ad — bc = 0. 
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Let’s write 


Since the determinant of A — uI is 0, we have ps — qr = 0. Now we define 


(=) 
x= ; 
=r 
Then 
E _{P 4 S\ (ps—aqr)\ _ (ps—ar 
(ASE 1) (i) E Gee E ( 0 | 
But we’ve just seen that ps — qr = 0, so 
(A — pI)x = 0. 


It follows that if x isn’t the zero vector, then it’s an eigenvector of A with 
corresponding eigenvalue p. 


If x is the zero vector, then we can define 


Ae 


and use a similar argument to see that if y isn’t the zero vector, then it’s 
an eigenvector of A with corresponding eigenvalue u. 


In the exceptional case when both x and y are equal to the zero vector, all 
four numbers p, q, r and s are 0, so 


0 0 
a-m- (i ve 


In this case, we have 
(A — pI)v =0 


for any vector v, so any non-zero vector is an eigenvector of A with 
corresponding eigenvalue p. 


To summarise, we saw earlier that any eigenvalue of A is a solution of the 
characteristic equation of A, and now we’ve seen that any solution of the 
characteristic equation of A is an eigenvalue of A. These important 
observations are encapsulated in the box below. 


Characteristic equation of a 2 x 2 matrix 


The characteristic equation of a 2 x 2 matrix A = ( J is the 
quadratic equation in À given by 
Gea ee ad te) 


The eigenvalues of A are the solutions of this equation. 


1 Eigenvalues and eigenvectors of matrices 


Soe : . {2 . 
For example, the characteristic equation of the matrix G s) is 


A? —5A+6=0. 


As you have seen, the eigenvalues — which are the solutions of the 
characteristic equation — are 2 and 3. 


In the next example and activity you’ll practise finding characteristic 
equations. Later on, you'll not only find characteristic equations, but also 
solve them to give the eigenvalues. 


Activity 4 Finding the characteristic equations of 2 x 2 matrices 


Find the characteristic equation of each of the following matrices. 
—2 1 -3 1 
(a) = i) (b) ( 4 k 


You may have noticed that the constant term of the characteristic 
equation of a 2 x 2 matrix is equal to the determinant of the matrix. That 
is, the matrix 


a b 
saa 
has characteristic equation 


à? — (a+ d)À +ad — bc = 0, 


and the constant term is the number ad — bc, which is the determinant 
of A. 
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The coefficient of À in the characteristic equation is —(a +d). You'll recall 
that the elements a and d are said to lie on the leading diagonal of the 
matrix (the diagonal that starts at the top-left element and ends at the 
bottom-right element). Their sum a + d is known as the trace of the 
matrix, so the coefficient of A in the characteristic equation of A is equal 
to minus the trace of A. 


Trace of a 2 x 2 matrix 


a bd 
Let a = (7 d 


It is the sum of the elements on the leading diagonal of A. 


) . The trace of A is the number a + d, written tr A. 


With the determinant and trace, there is now an alternative, equivalent 
form of the characteristic equation, which you may find easier to remember. 


Equivalent form of the characteristic equation 


The characteristic equation of a 2 x 2 matrix A is the quadratic 
equation in A given by 
V = (AA tA =A 


Example 3 Using the trace and determinant to find the 
characteristic equation of a 2 x 2 matrix 


Find the trace and determinant of the matrix 


1 -—2 
a 
and hence find its characteristic equation. 


Solution 
®. The trace and determinant are 

trA=a+d and detA=ad— be. 
e 
We have 

trA=1+5=6, 

det A = 1 x 5 — (—2) x 2 = 9. 

@. The characteristic equation is A? — (tr A)A + det A = 0. @ 


So the characteristic equation of A is 
M — 6A +9=0. 
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Activity 5 Using the trace and determinant to find the characteristic 
equations of 2 x 2 matrices 


For each of the following matrices, find the trace and determinant and 
hence find the characteristic equation. 


wa(h 4) wal) 


Once you’ve found the characteristic equation of a matrix, you can then 
find the eigenvalues of the matrix by solving the characteristic equation. 
This is a quadratic equation, so it either has one repeated solution or two 
distinct solutions. Each of these solutions is a complex number that may 
or may not be a real number. Since the solutions are the eigenvalues of the 
matrix, we have the following observation about the number of eigenvalues 
of a 2 x 2 matrix. 


Number of eigenvalues of a 2 x 2 matrix 


Every 2 x 2 matrix has either one or two eigenvalues, each of which is 
a complex number that may or may not be a real number. 


If the matrix has only one eigenvalue, then because it is a repeated 
solution of the characteristic equation, it is said to be 
a repeated eigenvalue of the matrix. 


For example, consider the characteristic equation A? — 6A + 9 = 0 that was 
obtained in Example 3. The quadratic expression can be factorised to give 


A—8/ =0. 


The number 3 is a repeated solution of this equation, so 3 is a repeated 
eigenvalue of the matrix. 


Let’s now practise finding eigenvalues of 2 x 2 matrices. The strategy 
you’ve learned for doing so is summarised in the next box. 


Strategy: 
To find the eigenvalues of a 2 x 2 matrix 


Let A = k a 

1. Find the trace and determinant of A, 
trA=a+d and detA=ad-— bc. 

2. Write down the characteristic equation 
d* — Gr A)A + det A = 0. 


3. Solve the characteristic equation to obtain the eigenvalues of A. 
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For example, you saw in Example 2 that the matrix 


(a) 


has characteristic equation A? — 4A +3 = 0. You can solve this equation by 
factorising the quadratic expression to give 


(A—1)(A— 3) =0. 
The eigenvalues are the solutions of this equation, namely 1 and 3. 


There is a useful way to check whether the values you obtain for the 
eigenvalues are correct. To understand how this check works, multiply out 
the expression (A — 1)(A — 3) on the left-hand side of the equation above to 
obtain 


Pasay). 


As you can see, the coefficient of À in this quadratic expression is equal to 
minus the sum of the eigenvalues. You learned earlier that the coefficient 
of À is also equal to minus the trace of the matrix. Therefore the sum of 
the eigenvalues is equal to the trace. 


You can use this observation to check that the eigenvalues you’ve found 
are correct: add them together and check that you obtain the trace. 


A useful way to check your eigenvalues 


The sum of the eigenvalues of a 2 x 2 matrix is equal to the trace of 
the matrix. 


If the matrix has only a single, repeated eigenvalue, then it must be 
counted twice in the sum. 


For instance, you saw in Example 3 that the matrix 


(5) 


has characteristic equation A? — 6\ + 9 = 0, and earlier we solved this 
equation to see that the matrix has a single, repeated eigenvalue 3. Since it 
is repeated, we count it twice when we work out the sum of the 
eigenvalues, to give 


34+3=6, 
which is equal to the trace of the matrix. 


This check doesn’t guarantee that the eigenvalues you’ve found are correct; 
after all, it may be that you’ve found two incorrect numbers for the 
eigenvalues that happen to add up to the trace. However, it’s still a useful 
quick way of testing for numerical mistakes. 
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There is another way to check whether the values you obtain for the 
eigenvalues are correct: their product must equal the determinant of the 
matrix. For example, the eigenvalues of the matrix 


(a) 


from Example 2 are 1 and 3, and their product is 3, which is equal to the 
determinant of the matrix. This check is less useful than the trace check, 
because it’s often more difficult to work out a product than a sum, and 
more difficult to work out a determinant than a trace. In this unit we use 
the trace test, but not the determinant test. 


Example 4 Finding the eigenvalues of 2 x 2 matrices 
Find the eigenvalues of the following matrices. 
(a) A= te J b AS E i) 
Solution 
(a) @ Find the trace and determinant of the matrix. © 
We have 
trA=14+2=3, 
det A = 1 x 2 — 2 x 3 = —4. 


@. Write down the characteristic equation 
d? — (tr A)A+ det A = 0. A 


So the characteristic equation of A is 
A -—3\-4=0. 


@. Solve the characteristic equation to find the eigenvalues. This 
one can be solved by factorising the quadratic expression. © 


Hence 
(A+ 1)(A— 4) =0, 
so the eigenvalues of A are —1 and 4. 
@. Check that the sum of the eigenvalues is equal to the trace. ® 
(Check: (—1) +4=3=trA.) 
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(b) &. Find the trace and determinant of the matrix. © 
We have 
cA 837l = 4 
det A = 3 x 1 — 1 x (—2) = 5. 


@®. Write down the characteristic equation 
7 — (trA)A\+detA =0. # 


So the characteristic equation of A is 
PFs eal) 


@. Solve the characteristic equation to find the eigenvalues. Use 
the quadratic formula. ©& 


Hence 
n —(—4) + y (—4) —4x1x5 
E Dy il 
447-4 
E 2 
d% 
se 
= 2 se 4. 


So the eigenvalues of A are 2 — i and 2 + i. 
@®. Check that the sum of the eigenvalues is equal to the trace. & 
(Check: (2 — i) + (2+74) =4=trA.) 


Activity © Finding the eigenvalues of 2 x 2 matrices 


Find the eigenvalues of the following matrices. 


@ A=(T 4) Ma=(3 4) ©4-( 4) 


1.3 Finding eigenvectors 


After you’ve found the eigenvalues of a 2 x 2 matrix by solving the 
characteristic equation, you can then find the corresponding eigenvectors 
by using the eigenvector equations. In this subsection you’ll see how to do 
this for real eigenvalues (eigenvalues that are real numbers). The 
strategy for finding eigenvectors is similar when the eigenvalues are not 
real numbers. 


112 


1 Eigenvalues and eigenvectors of matrices 


Let’s demonstrate the strategy using the matrix 


w 


In Example 4(a) you learned that A has two eigenvalues, —1 and 4. Near 
the start of the previous subsection you saw that the eigenvector 

equation Ax = Ax can be stated in the equivalent form (A — AI)x = 0. In 
the example here, the equivalent form of the eigenvector equation is 


1—à 2 z\ /0 
3 2— AJ \y}  \07? 
where A is an eigenvalue of A, and x = C) 


When à = —1, we obtain 


2 2\/x\ _ (0 
aay os RO 
We can write this matrix equation as a pair of simultaneous equations, 


2x + 2y = 0, 
3x + 3y = 0. 
By dividing both sides of the first equation by 2, and dividing both sides of 


the second equation by 3, you see that each of the two equations is 
equivalent to the single equation 


xty=0; thatis, y=-—z. 
This final equation y = —2 is the equation of the line shown in Figure 5. 


The position vector x of any point on this line satisfies (A — AI)x = 0, 
with À = —1, so any such vector other than the zero vector is an 
eigenvector of A corresponding to the eigenvalue —1. For example, 
choosing x = 1 on the line y = —x gives y = —1, so 


(i) 


is an eigenvector of A corresponding to the eigenvalue —1. Similarly, 
choosing x equal to 5 and —2 in turn on the line y = —2 gives vectors 


(3) (4) 


each an eigenvector of A corresponding to the eigenvalue —1. More 
generally, any eigenvector of A corresponding to the eigenvalue —1 has the 
form 


(3 


for some non-zero number z. 


Sy 


y=-r 


Figure 5 The line y = =x 
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You can check that such vectors are indeed eigenvectors of A corresponding 
to the eigenvalue —1 by verifying that the eigenvector equation 


G 2) (a)=-G) 


is satisfied. Geometrically, this equation demonstrates that each point on 


the line y = —< is scaled by the linear transformation represented by A 
by the factor —1, so the line y = —~ is an invariant line of the 
transformation. 


This method for finding eigenvectors is summarised in the following 
strategy box. 


Strategy: 
To find an eigenvector of a 2 x 2 matrix corresponding to a 
given eigenvalue 


To find an eigenvector of a 2 x 2 matrix A corresponding to an 
eigenvalue A of A, carry out the following steps. 


1. Write down the eigenvector equation (A — AI)x = 0, 


where x = C) is an eigenvector to be determined. 


2. Rewrite this equation as a pair of simultaneous equations 
in x and y. 


3. These simultaneous equations are equivalent to a single equation 
in z and y. Any solution of this equation other than z = y = 0 
gives an eigenvector of A corresponding to the eigenvalue A. 
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@. Write this matrix equation as a pair of simultaneous equations. © 
This gives 
—3z + 2y = 0, 
Sp — 7) = (0). 
@. Find a solution of these equations. ® 
This pair of equations is equivalent to the single equation 
2) = See. 


@®. Choose a value of x that gives a simple solution, preferably 
without fractions. Choosing x = 1 would give y = 3, which involves a 
fraction. It’s better to choose x = 2 because then y = 3, and neither x 
nor y is a fraction. @ 


If x = 2, then 2y = 3 x 2, so y= 3. Hence (a is an eigenvector of 


the matrix corresponding to the eigenvalue 4. 


@. Check your answer by making sure that the eigenvector 
equation Ax = Ax is satisfied. ® 


ome (JOO 


Activity 7 Finding eigenvectors of a 2 x 2 matrix 


The matrix ¢ >) has eigenvalues —2 and 5. For each of these 


eigenvalues, find a corresponding eigenvector. 


As you’ve seen, to find an eigenvector of a 2 x 2 matrix A corresponding 
to an eigenvalue À of A, you must find a solution other than x = y = 0 of a 
pair of simultaneous equations 

pz + qy = 0, 

rx + sy =0, 
where p, q, r and s are real numbers (they are the elements of the 
matrix A — AI). For instance, in Example 5 this pair of simultaneous 
equations is 

—3x2 + 2y = 0, 

3x — 2y = 0. 
If p and q are not both 0, then the equation px + qy = 0 represents a line 
in the plane that passes through the origin. If p and q are both 0, then the 
equation is 

Ox + Oy = 0. 
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Figure 6 The line y = žr 
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This equation is satisfied by every pair of values of x and y, so it 
represents the entire plane. Of course, similar comments apply to the 
line rx + sy = 0. 


Often when calculating eigenvectors by solving the pair of simultaneous 
equations px + qy = 0 and rx + sy = 0, you'll find that neither equation 

is Ox + Oy = 0, so both equations represent lines through the origin. In 
fact, they must represent the same line through the origin because 
otherwise the two lines would intersect only at the origin, which would 
imply that the only solution of the simultaneous equations is x = y = 0. If 
you discover that the only solution is x = y = 0, then you’ve made a 
mistake, because you should be able to find a (non-zero) eigenvector from 
one of the solutions. 


As an example of when both equations represent the same line, consider 
again the simultaneous equations 

—3x + 2y = 0, 

3x — 2y = 0, 
from Example 5. Each of these equations is equivalent to the single 
equation 2y = 3a; that is, y = 3x. The line given by this equation is shown 
in Figure 6. 
Sometimes you'll find that one of the two equations px + qy = 0 
and ra + sy = 0 represents a line, but the coefficients of the other equation 
are both 0, so it represents the entire plane. For example, the pair of 
simultaneous equations 

2r +y=0, 

Ox + Oy = 0, 
is of this type. In this case, the solutions of the simultaneous equations are 


given by the equation representing a line, because all values of x and y 
satisfy the equation with 0 coefficients. 


There is one final very special pair of simultaneous equations that we 
haven’t yet considered, when all the coefficients of the simultaneous 
equations are 0. To understand this case, let’s suppose that we are trying 
to find an eigenvector of a 2 x 2 matrix A corresponding to an eigenvalue A 
of A. In this very special case, when we write the eigenvector equation 


(A — XI)x =0 


as a pair of simultaneous equations, all the coefficients are 0. It follows 
that A — AI must be the zero matrix; that is 


0 0 
awe (2%). 


Adding AI to both sides of this equation gives 
A=AI. 


1 Eigenvalues and eigenvectors of matrices 


So A is a scalar multiple of the identity matrix. As you saw in Unit 6, 
matrices of this type represent linear transformations called dilations. For 
example, the matrix 


sGD 


represents a dilation. This matrix has an eigenvalue 5 because the 
eigenvector equation 


Ax = 5x 


is satisfied for any vector x whatsoever. What is more, because this 
equation is satisfied for every vector x, the number 5 is the only eigenvalue 
of A. 


This discussion about the eigenvectors corresponding to a real eigenvalue is 
summarised below. 


Eigenvectors corresponding to a real eigenvalue 
Suppose that A is a 2 x 2 matrix with a real eigenvalue A. 


e If A ATI, then the eigenvectors of A corresponding to À are the 
position vectors of the points on a single line through the origin, 
other than the position vector of the origin itself. 


e If A = AI, then every non-zero vector is an eigenvector of A 
corresponding to X. 


Let’s now see an example of finding eigenvalues and eigenvectors of a 
matrix. 
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@. Write down the characteristic equation 
—(rA)\ det A=0. # 


So the characteristic equation of A is 
MW =-444=0. 


®. Solve the characteristic equation to find the eigenvalues. This one 
can be solved by factorising the quadratic expression. © 


Hence 
(A z 2)? z= 0, 
so A has only a single, repeated eigenvalue 2. 


@. You should now carry out a quick mental check that the sum of 
the eigenvalues is equal to the trace of A. Later on, after you’ve 
found an eigenvector of A corresponding to the eigenvalue 2, you’ll 
check both this eigenvalue and eigenvector using the eigenvector 
equation for 2. ® 


@. Find the eigenvectors. To do this, first write down the eigenvector 
equations of the matrix in the form (A — \I)x = 0. ® 


The eigenvector equations have the form 


Ua" 523) G) =): 


Substituting A = 2 gives 


(=) G) =): 


@. Write this matrix equation as a pair of simultaneous equations. © 


This gives 
T Y= 0, 
—x#+y=0. 


@. Find a solution of these equations. ® 
This pair of equations is equivalent to the single equation 
h = i- 


@. Choose a value of x that gives a simple solution. ® 


ie =i then y = ii Hence G 


i) is an eigenvector of A corresponding 


to the eigenvalue 2. 


@. Check your answer by making sure that the eigenvector 
equation Ax = Ax is satisfied. ® 


a(i 0) -@)-20) 
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Activity 8 Finding eigenvalues and eigenvectors of 2 x 2 matrices 


Find the eigenvalues of the following matrices, and for each eigenvalue find 
a corresponding eigenvector. 


(a) A= t >) (b) A= € s) (c) A= & 2 
@ a=(4 o 


All the matrices that you’ll meet in this unit are real matrices, which 
means that their elements are real numbers. As you’ve seen already, a real 
matrix may have eigenvalues that are not real numbers. For example, you 
saw in Example 4(b) that the matrix 


(2) 


has eigenvalues 2 — 7 and 2+ i. There are eigenvectors corresponding to 
these eigenvalues, but the components of those eigenvectors are complex 
numbers that are not real numbers. Here we restrict our attention to real 
eigenvectors (eigenvectors whose components are real numbers) and do 
not consider eigenvectors corresponding to eigenvalues that are not real 
numbers. 


1.4 Matrices larger than 2 x 2 


Many of the techniques that you’ve learned for finding eigenvalues and 
eigenvectors of 2 x 2 matrices can also be applied to bigger square 
matrices. In particular, the method used to find the characteristic equation 
of a 2 x 2 matrix can be adapted for n x n matrices. The result of 
applying this method to an n x n matrix A is a polynomial equation, 
again called the characteristic equation of A, whose solutions are the 
eigenvalues of A. 


The degree of the characteristic equation of an n x n matrix is n. A 
polynomial equation of degree n has at least 1 and at most n solutions. It 
follows that an n x n matrix has at least 1 and at most n eigenvalues. 


For example, the characteristic equation of a 2 x 2 matrix is a quadratic 
equation, as you’ve seen already, so a 2 x 2 matrix has 1 or 2 eigenvalues. 
As another example, consider the 3 x 3 matrix 


0 1 1 
A=j]|-1l 2 1 
—1 -1 4 
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Using the computer algebra system (see Activity 10) you can show that 
the characteristic equation of A is the cubic equation 


dM + 6d? — 111 +6 =0. 


In general, cubic equations are difficult to solve; however, this particular 
equation happens to have exactly three solutions, namely 1, 2 and 3 
(which you can check yourself). These are the eigenvalues of A. 


You can find eigenvectors of A corresponding to the eigenvalues 1, 2 and 3 
by using similar methods to those that you used for 2 x 2 matrices. The 
difference this time is that you have to solve three simultaneous equations 
in three unknowns, rather than two simultaneous equations in two 
unknowns. Although you can do this by hand, the procedure takes a lot 
longer than it does for 2 x 2 matrices. 


In the next activity, you are asked to check that three vectors are 
eigenvectors of A, and to find the eigenvalues corresponding to each of 
them. 


Activity 9 Finding the eigenvalues corresponding to given 
eigenvectors 


Show that the vectors 


1 1 2 
1], 1 and 1 
1 2 1 
are eigenvectors of the matrix 
0 1 1 
—-l1 21], 
—1 -1 4 


and for each eigenvector find the corresponding eigenvalue. 


Usually the most practical way of finding the eigenvalues and eigenvectors 
of an n x n matrix, where n > 3, is to use the computer algebra system. 


TU Activity 10 Finding eigenvalues and eigenvectors with a computer 

Work through Section 8 of the Computer algebra guide, where you will 
learn how to find eigenvalues, eigenvectors and characteristic equations of 
matrices using the computer algebra system. 


The matrices in the remainder of this unit are all of size 2 x 2. However, 
many (but not all) of the techniques that you’ll meet can also be applied 
to bigger square matrices. 
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Olga Taussky-Todd: a torchbearer for matrix theory 


The mathematician Olga Taussky-Todd was born in 1906 in the city 
of Olomouc, which at the beginning of the twentieth century was part 
of the Austro-Hungarian Empire, but now belongs to the Czech 
Republic. She moved to England in 1935, and during the Second 
World War she worked at the National Physical Laboratory in 
Teddington, Middlesex. Part of her work there involved using 
matrices to analyse the stability of aircraft. 


One of her major mathematical achievements was to develop a result 
about eigenvalues known as the Gershgorin circle theorem, and utilise 
it in her work on aeronautics. Later in life she wrote that on first 
meeting the result, she ‘immediately tinkered with the theorem, 
applying it to a very nasty looking 6 x 6 matrix’. She published her 
findings in a report to the Aeronautical Research Council in 1947. 


Because of her contributions to the theory of matrices, the 
British-American mathematician Hans Schneider wrote that 
Tausky-Todd’s ‘work has altered the consciousness of several 
generations of matrix theorists’. 


(Source: Schneider, H. (1977) ‘Olga Taussky-Todd’s influence on 
matrix theory and matrix theorists’, Linear and Multilinear Algebra, 
vol. 5, no. 3, pp. 197-224) 


2 Eigenvalues and eigenvectors of 
special types of matrices 


Using the strategies of the previous section, you can find the eigenvalues 
and eigenvectors of any 2 x 2 matrix. However, for some types of matrices 
there are far quicker ways of identifying the eigenvalues and eigenvectors. 
Here you'll learn about some of these special types of matrices, and you'll 
also study the eigenvalues and eigenvectors of matrices representing many 
of the standard types of linear transformations that you met in Unit 6, 
such as flattenings, reflections, rotations, scalings and shears. 


2.1 Triangular matrices 


In this subsection you'll learn about a class of matrices called the 
triangular matrices. The eigenvalues of matrices from this class are 
particularly easy to identify. 


Olga Taussky-Todd (1906-95) 
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Triangular matrices 


A 2 x 2 triangular matrix is a matrix of one of the forms 
a b a 0 
Oy Pe wy 


Matrices of the first type are called upper triangular matrices and 
matrices of the second type are called lower triangular matrices. 


For example, 


(5 5)» (ra) om (0 9) 


are all triangular matrices. The first matrix is an upper triangular matrix, 
the second is a lower triangular matrix, and the third is both an upper and 
a lower triangular matrix. In contrast, 


0 2 d -1 7 
ga = 4 0 
are not triangular matrices. 


Notice that 2 x 2 diagonal matrices, that is, matrices of the form 


(0 4) 


are triangular matrices. 


The eigenvalues of a triangular matrix are easy to identify using the 
following observation. 


Eigenvalues of 2 x 2 triangular matrices 


The eigenvalues of a 2 x 2 triangular matrix are the elements on the 
leading diagonal. 


For example, the eigenvalues of the triangular matrix 


e 


are 2 and 9. To see why this is so, remember that the eigenvalues of A are 
the solutions of the characteristic equation of A, and the characteristic 
equation is obtained by writing the left-hand side of the equation 


det(A — AT) = 0 


as a quadratic expression in À. In this case 
2-r 3 
au ( 0 9- N ? 
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and the determinant of this matrix is (2 — A)(9 — A). Therefore the 
characteristic equation is 


G=- 


This equation has solutions 2 and 9, which are the elements on the leading 
diagonal of A. 


You can also quickly identify an eigenvector of a 2 x 2 triangular matrix, 
using the following observation. 


Eigenvectors of 2 x 2 triangular matrices 


If A= 6 i) , then (a) is an eigenvector of A corresponding to the 


eigenvalue a. 


If A= ( , then (7) is an eigenvector of A corresponding to the 


eigenvalue d. 


Let’s see why this is true for upper triangular matrices; you can reason in a 
similar way for lower triangular matrices. Observe that 


WORO RO] 


1\. : . . : 
so G is an eigenvector of the matrix corresponding to the eigenvalue a. 


To find an eigenvector of the matrix corresponding to the other 
eigenvalue d you should use the strategy given in the previous section. 


As a 2 x 2 diagonal matrix is both an upper and a lower triangular matrix, 


the vectors (o) and @ 


1) are both eigenvectors of such a matrix. 


Eigenvalues and eigenvectors of 2 x 2 diagonal matrices 


Let A = ( 7 . The eigenvalues of A are the elements a and d on 
the leading diagonal, and E and a are eigenvectors of A 


corresponding to the eigenvalues a and d, respectively. 
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Example 7 Finding eigenvalues and eigenvectors of 2 x 2 triangular 
matrices 

Find the eigenvalues of the following matrices, and for each eigenvalue 
find a corresponding eigenvector. 


@ (24) © (3) 


Solution 


(a) @. This matrix is a triangular matrix, so the eigenvalues are the 
elements on the leading diagonal. © 


Since the matrix is a triangular matrix, the eigenvalues are the 
elements on the leading diagonal, namely 1 and 4. 

@. Next find the eigenvectors. The matrix is a lower triangular 
: is an eigenvector corresponding to the 
eigenvalue given by the bottom-right element of the matrix. © 


matrix, so the vector 


The vector G) is an eigenvector of the matrix corresponding to 
the eigenvalue 4. 


@. Find an eigenvector corresponding to the other eigenvalue in 
the usual way, using the corresponding eigenvector equation. © 


The eigenvector equations have the form 
1-A 0 ai (0 
2 AAAA T Oo 


When à = 1, we obtain 


e)0-0 


@. Write this matrix equation as a pair of simultaneous 
equations. ® 


This gives 
Oz + Oy = 0, 
22 + 3y = 0. 


This pair of equations is equivalent to the single equation 
2x+3y=0; thatis, 3y = —2z. 


@. Choose a value of x that gives a simple solution. ® 


If x = 3, then 3y = —2 x 3, so y = —2. Hence @ is an 


eigenvector of the matrix corresponding to the eigenvalue 1. 


2 Eigenvalues and eigenvectors of special types of matrices 


Activity 11 Finding eigenvalues and eigenvectors of 2 x 2 triangular 
matrices 


Find the eigenvalues of the following matrices, and for each eigenvalue find 
a corresponding eigenvector. 


oC) aCA of} of 2) 


In Unit 6 you learned about a class of linear transformations called 
shears. In particular, you saw that horizontal shears and vertical 
shears are transformations represented by triangular matrices of the forms 


1k 1 0 
(0 1) amt (& 4) 


respectively. 


Activity 12 Finding eigenvalues and eigenvectors of a matrix 
representing a horizontal shear 


Find the eigenvalues of the matrix 


(0 i): 


and for each eigenvalue find a corresponding eigenvector. 
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Suppose that A is the matrix 


(0 3) 


from the preceding activity, and let f be the horizontal shear represented 
by A. Let i and j be the usual unit vectors (a) and Gr and let i’ = Ai 
and j/ = Aj. The effect of f on the unit square is shown in Figure 7. 


YA YA 


Ry 
E 
—s, 
= 
=] 
M 

Sv 


Figure 7 ‘The effect of a horizontal shear 


You learned in Unit 6 that a horizontal shear fixes every point on 

the x-axis, and moves points that are not on the x-axis horizontally. In 
particular, the shear f fixes the point (1,0), as you can see from Figure 7. 
Using matrix notation, we can write the observation that f(1,0) = (1,0) as 


(a) = (0) 


This gives us a geometric way of seeing that is an eigenvector of A 


1 
0 
corresponding to the eigenvalue 1, a fact that you may already have 
verified in Activity 12. 


However, the matrix A has no eigenvectors other than scalar multiples 
0 


horizontally by f, so the image point cannot be a scalar multiple of the 
original point. 


of a): This is because any point that doesn’t lie on the x-axis is moved 


Similar remarks about eigenvalues and eigenvectors can be made for other 
horizontal and vertical shears. We summarise them below. 


2 Eigenvalues and eigenvectors of special types of matrices 


Eigenvalues and eigenvectors of matrices representing 
horizontal and vertical shears 


The matrices 


i & i @ 
€ ') and G a where k £0, 


which represent a horizontal shear and a vertical shear, respectively, 
each have only a single, repeated eigenvalue 1. The corresponding 
eigenvectors are all non-zero scalar multiples of 


(c) = G) 


respectively. 


2.2 Flattenings 


In Unit 6 you learned that a linear transformation is a flattening if its 
image set is either a line through the origin, or, in the case of the zero 
transformation, the set containing the origin alone. You also saw that to 
test whether a linear transformation is a flattening you can check whether 
the determinant of the matrix representing the transformation is zero. 
Here we investigate the eigenvalues and eigenvectors of matrices 
representing flattenings. 


Suppose that the 2 x 2 matrix A represents a flattening, so det A = 0. The 
characteristic equation of A then becomes 


X? — (tr A= 0 
The quadratic expression on the left-hand side factorises to give 
A(A — tr A) = 0, 


so the eigenvalues of A are 0 and tr A. 


Eigenvalues of 2 x 2 matrices with determinant zero 


The eigenvalues of a 2 x 2 matrix A with determinant zero 
are 0 and trA. 


If tr A = 0, then A has only a single, repeated eigenvalue 0. 


For example, the matrix 


0) 


has determinant zero, and tr A = 5, so its eigenvalues are 0 and 5. 
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Figure 8 The 
transformation f flattens the 
plane onto the line y = 3a 
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As the determinant of this matrix A is zero, the linear transformation f 
represented by A is a flattening. It’s not the zero transformation, so the 
image set of f is a line through the origin. We can find out what that line 
is by looking at the images of points under f. 


For example, under f the images of the points with position vectors 


=() m0) 


must both be points on the line. But Ai and Aj are just the first and 
second columns of A, respectively, namely 


(3) ea (a) 


Each of these vectors is a scalar multiple of , so we see that f flattens 


1 
3 
the plane onto the line that passes through the origin and the point (1,3). 
This is the line y = 3a shown in Figure 8. 


The line y = 3x is mapped to itself by f, so it’s an invariant line of f. It 
follows that any point on this line is mapped by f to a scalar multiple of 
itself, so the position vector of any non-zero point on the line is an 
eigenvector of A. In particular, the two columns of A are eigenvectors 
of A. Let’s check this for the second column: 


a(3)= (6 3) G) = Gas) 563) 


Therefore t 


(and so is (6) 


Notice that the eigenvectors given by the columns of A both correspond to 
the eigenvalue tr A = 5, rather than to the other eigenvalue 0. This is 
because they are the position vectors of points on the line y = 3x, which is 
mapped to itself by f, rather than to the origin. 


) is an eigenvector of A corresponding to the eigenvalue 5 


With similar reasoning we can obtain the following observations about the 
eigenvectors of a 2 x 2 matrix with determinant zero. 


Eigenvectors of 2 x 2 matrices with determinant zero 


Let A be a 2 x 2 matrix with determinant zero. Each column of A, if 
non-zero, is an eigenvector of A corresponding to the eigenvalue tr A. 


For example, the determinant and trace of the matrix 


(i =i) 


are both zero, so it has only a single, repeated eigenvalue 0. 


2 Eigenvalues and eigenvectors of special types of matrices 


Each column of the matrix is 


which is an eigenvector of the matrix corresponding to the eigenvalue 0. 


In this example, the observations you’ve seen about eigenvalues and 
eigenvectors of 2 x 2 matrices with determinant zero allow you to quickly 
find the eigenvalue and a corresponding eigenvector of the matrix. This 
matrix is special though, because both the determinant and the trace 
are zero, so there is only one eigenvalue, namely 0. 


A 2x 2 matrix A whose determinant is zero but trace is not zero has two 
eigenvalues, 0 and tr A, and the observations that you’ve met so far in this 
subsection don’t help you to find an eigenvector of A corresponding to the 
eigenvalue 0. So in many cases, the observations won’t save you much 
time, and rather than memorising them you may prefer to stick with the 
general strategies for finding eigenvalues and eigenvectors given in the 
previous section (and use these observations as additional checks). 


Example 8 Finding eigenvalues and eigenvectors of a 2 x 2 matrix 
with determinant zero 
Find the eigenvalues of the matrix 


4 6 
a=( a) 
and for each eigenvalue find a corresponding eigenvector. 


Solution 
@. To find the eigenvalues of the matrix, first find the trace and 
determinant. ® 
We have 
trA =4+3=7, 

det A =4x 3-6x2=0. 
@. Since the determinant is zero, the eigenvalues are equal to 0 
and tr A. You don’t have to remember this observation though; it 


doesn’t take long to go through the usual strategy of finding 
eigenvalues by solving the characteristic equation. © 


So the characteristic equation of A is 


v= 7N=0: 
Hence 
MA — 7) = 0, 


so the eigenvalues of A are 0 and 7. 


129 


Unit 11 Eigenvalues and eigenvectors 


130 


@. Next find the eigenvectors. Write down the eigenvector equations 
of the matrix in the form (A — \I)x = 0. & 


The eigenvector equations have the form 


4— 6 He o fa) 
2° B= Ni yp Oy 
@. Substitute A = 0 and solve the resulting equations to find an 


eigenvector corresponding to the eigenvalue 0. © 


When A = 0, we obtain 


(2 3) ()=(0). 


This gives 
4x + 6y = 0, 
P => dy) = 0. 


This pair of equations is equivalent to the single equation 


27+ 3y=0; thatis, 3y = —2z. 


If x = 3, then 3y = —2 x 3, so y = —2. Hence o is an eigenvector 


of A corresponding to the eigenvalue 0. 


a(i D-O- 


@. Solve the eigenvector equation with À = 7 to find an eigenvector 
of A corresponding to the eigenvalue 7. © 


When A = 7, the eigenvector equation is 
=o 2 OV (av 0 
2 =—4) \y)p ~ \O) 
This gives 


—3x + ôy = 0, 
2x — 4y = 0. 


This pair of equations is equivalent to the single equation 
m= =0 a is, a = Zo. 
2 5 
If y = 1, then z = 2. Hence Gl is an eigenvector of A corresponding 


to the eigenvalue 7. 


come ($9)()-(3)-10) 


2 Eigenvalues and eigenvectors of special types of matrices 


The matrix 


a(S) 


from the preceding example represents a flattening f onto a line that passes 
through the origin. Since the columns of A are position vectors of points 
in the image of f, and both columns are scalar multiples of the eigenvector 


oj 


we see that the image line passes through the point (2,1). This line has 
equation y = iz. 


3 f 
_9 corresponding 
to the eigenvalue 0. Since non-zero scalar multiples of this eigenvector are 
also eigenvectors of A corresponding to the eigenvalue 0, we see that f 
maps the whole of the line y = -2r that passes through the origin and the 
point (3, —2) to the origin, as shown in Figure 9. 


The other eigenvector we found in the example was 


image of y = -ĉr 


Figure 9 The transformation f flattens the plane onto the line y = ix 
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Figure 10 Rotation through 
an angle of 7/2 about the 


origin 
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Activity 13 Finding eigenvalues and eigenvectors of 2 x 2 matrices 
with determinant zero 


Find the eigenvalues of the following matrices, and for each eigenvalue find 
a corresponding eigenvector. 


oa(1 i) ma) 


2.3 Rotations and reflections 


In this subsection, you’ll study the eigenvalues and eigenvectors 

of 2 x 2 matrices that represent rotations and reflections. You’ll see that 
matrices representing rotations about the origin typically have no real 
eigenvalues, whereas matrices representing reflections in lines through the 
origin each have two real eigenvalues. 


Rotations 


In Unit 6 you learned that a rotation through an angle 0 about the origin 
is represented by the matrix 


cos@ —sind 

sin 6 cos 6} ` 
For example, the rotation through an angle 7/2 about the origin is 
represented by the matrix 


0 —1 

1 0J’ 
This rotation is illustrated in Figure 10, which shows the image of the unit 
square after applying the rotation. 


In the next activity you'll see that this matrix has no real eigenvalues. 


Activity 14 Finding the eigenvalues of a matrix representing a 
rotation 


Find the eigenvalues of the matrix 
0 -1 
te a o) 


The matrix representing a rotation from the preceding activity has no real 
eigenvalues, so it also has no real eigenvectors. This is to be expected, for 
if there was a real eigenvector, then it would be the position vector of a 
point that lies on an invariant line of the rotation that passes through the 
origin. But there is no such invariant line for a rotation by 7/2. 


2 Eigenvalues and eigenvectors of special types of matrices 


Reasoning in a similar way, you can see that almost any matrix that 
represents a rotation about the origin has no real eigenvalues. There are 
two exceptions to this, namely the matrices 


1 0 -1 0 
t= (4 i and -1=( 0 =) 


which represent rotations about the origin through angles 0 and 7, 
respectively. As you can see, both matrices are scalar multiples of the 
identity matrix. The matrix I has only a single, repeated eigenvalue 1, 
and —I has only a single, repeated eigenvalue —1, and every non-zero 
vector is an eigenvector of both matrices. 


Eigenvalues and eigenvectors of 2 x 2 matrices representing 
rotations 


With the exception of the matrices I and —I, every 2 x 2 matrix that 
represents a rotation about the origin has no real eigenvalues or real 
eigenvectors. 


Reflections 


In Unit 6 you saw that a reflection in a line @ through the origin is 
represented by the matrix 


e 2a sin 
A=|{. ; 
sin2a —cos2a 


where a is the angle of inclination of £. Figure 11 shows the effect of this 
reflection on the unit square. 


The matrices of such reflections, unlike those of rotations about the origin, 
always have real eigenvalues and eigenvectors. To see why this is so, notice 
that any point on the line @ is fixed by the reflection, so the position vector 
of that point must be an eigenvector of A with eigenvalue 1. For instance, 
the point (cos q, sin a) lies on £Z, as you can see from Figure 12, so we 
should expect that 


A J E o) l 
sina sin @ 
To verify this, we need to recall some trigonometric identities. 
Angle difference identities 


sin(A — B) = sin Acos B — cos Asin B 
cos(A — B) = cos Acos B + sin Asin B 


Sv 


Figure 11 Reflection in the 
line £ through the origin 


YA 


(cos a, sin a) 


COS Q x 


Figure 12 The 
point (cos qa, sinq) lies on £ 
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Figure 13 
reflection 
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Eigenvectors of a 


Now 


A (C82) _ (008 2a cosa + sin 2a sin a 
sin @ sin 2a cos a — cos 2a sina j ’ 


and using the angle difference identities this vector becomes 
cos(2a—a)\ _ (cosa 
sin(2Qa—a)})  \sina/’ 

as expected. 


Next consider a vector that is perpendicular to fe c) such 


cos @ 
their scalar product is 0. This new vector is the position vector of a point 
on a line through the origin that is perpendicular to ¢. After a reflection 
in £, such a point is scaled by the factor —1, so we should expect that 


— sina —sina 
A =— ; 
cos @ cos @ 


Again, you can check that this is true using the angle difference identities. 


as @ si") You know that these two vectors are perpendicular because 


The two eigenvectors of A that we’ve found are shown in Figure 13. 


Eigenvalues and eigenvectors of 2 x 2 matrices representing 
reflections 


The matrix 
cos 2q sin 2q 
sin2a —cos2a}’ 
which represents a reflection in the line through the origin with angle 


of inclination a, has eigenvalues —1 and 1, and corresponding 
eigenvectors 


— sina Cos Q 
and : ; 
COS @ sin @ 


respectively. 


In part (a) of the next activity you’ll revisit a matrix representing a 
reflection that you met right at the start of this unit. 


Activity 15 Finding eigenvalues and eigenvectors of matrices 
representing reflections 


Find the eigenvalues of the following matrices, and for each eigenvalue find 
a corresponding eigenvector. 


TORE 
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2.4 Generalised scalings 


Let’s begin by revisiting another matrix that you met at the start of this 
unit, namely 


azh n 


This represents a (2,3)-scaling, which scales horizontal distances by the 
factor 2 and vertical distances by the factor 3. You saw that 


(a) is an eigenvector of A corresponding to the eigenvalue 2, 


and 


(1) is an eigenvector of A corresponding to the eigenvalue 3. 


These two eigenvectors are perpendicular to one another. Here you’ll meet 
a more general class of scalings, in which each scaling is represented by a 
matrix that has two eigenvectors that, although not necessarily 
perpendicular, are not scalar multiples of one another. 


To introduce you to this more general class of scalings, we'll consider the 
matrix 


0 2 
ne(22) 
Since 
2 2 1 1 
aG)=G) = 4G)-20). 
you can see that A has eigenvalues 1 and 2, with corresponding 


eigenvectors (7) and Gi respectively. Let f be the linear 


transformation represented by A. The effect of f on the two eigenvectors 
is shown in Figure 14. 


YA y 


Ry 
Ry 


Figure 14 The effect of f on two eigenvectors 
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YA 
(1,1 
(2,1) 
(0, 0) 5 
Figure 15 Parallelogram 


with vertices (0,0), (2,1), 
(3,2) and (1, 1) 
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As you know, any point on the line that passes through the origin and the 
point (1,1) is scaled under f by the factor 2, because the position vector of 
any point on that line (other than the origin) is an eigenvector of A 

corresponding to the eigenvalue 2. With similar reasoning, you can see that 
every point on the line that passes through the origin and the point (2,1) is 
scaled under f by the factor 1; in other words, any such point is fixed by f. 


Let’s now try to work out how f transforms points that don’t lie on one of 
these two lines. To do this, consider the parallelogram with vertices 


(0,0), (2,1), (3,2) and (1,1). 


The position vectors of the second and fourth of these vertices are the two 
eigenvectors that we found. The parallelogram is shaded blue in Figure 15. 


Now ‘tile’ the rest of the plane using this base parallelogram to give a grid 
of parallelograms, as shown in Figure 16. Each line of the grid is parallel 
to one of the sides of the parallelogram. 


YA 


Figure 16 Grid of parallelograms 


Consider a vertex of this grid, such as the point (7,5) shown in Figure 16. 
You can write the position vector of this vertex as a sum of integer 


multiples of the eigenvectors ©) and ia). In this case, 
7\ _ (4 4 3 
5) 2 3 
2 1 
=2 ©) +3 C) 


2 Eigenvalues and eigenvectors of special types of matrices 


; ; ‘oe 
We can use this expression for (5) in terms of the eigenvectors of A to 


7 A f . 
determine A (5). as follows. First, using standard properties of matrices, 


we have 


(5) =4(2(2) +9(2)) =24(2) a 


But we already know that 
2 2 1 1 
a()=() wm aG)=2(), 
7 2 1 
a (3) a +3x2 (1): 


Simplifying the right-hand side gives 


ala) = C3): 


So f maps the point (7,5) to the point (10,8). 


SO 


In summary, to find the image of the point (7,5) under f you write 


í 2 1 
()=26) +9): 
and then observe that 
ih 2 1 
a(s)=2() +6), 
The first of these equations shows how to write (7) as the sum of 


a vector with direction (7) and a vector with direction (7) i 


The second of these equations shows that the image vector A (7) is 
i . i T 2 
obtained by scaling the part of the sum with direction (7) by the factor 1, 


and scaling the part with direction C) by the factor 2. 


You can find the image of any vertex of the grid under f in this way, by 
writing the position vector of that vertex as a sum of integer multiples of 
the eigenvectors 


(i) e G) 


in that order, and then scaling the first part of the sum by the factor 1 and 
scaling the second part of the sum by the factor 2. In fact, you can find 
the image of any point in the plane under f by writing the position vector 
of that point as a sum of real number (not necessarily integer) multiples of 
the two eigenvectors, and then scaling each part separately. 
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‘Obviously the architect did 
some tweaking but the basic 
idea was mine.’ 
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Try this strategy for working out the images of points in the next activity. 


Activity 16 Finding images of points under a linear transformation 


Find the images of the following points under the linear transformation f 
from above, which is represented by the matrix 


-o 


(a) (4,4) (b) (—2,-1) wW (1,—2) ik e =3 (7) =g C) J 


Figure 17 shows the effect of the transformation 


x 0 2 x 
(= Ci dG) 
that we’ve been discussing on the base parallelogram and the grid of 


parallelograms. 


YA YA 


Ry 
8 


Figure 17 The effect of a generalised scaling 


This transformation f is called a generalised scaling because of the way 
you can work out images of points by scaling in the directions of the two 
eigenvectors. In fact, you can find images of points in a similar way for any 
linear transformation of the plane that is represented by a matrix that has 
two real eigenvectors that are not scalar multiples of one another, so all 
such transformations are referred to as generalised scalings. 


Generalised scalings 


A generalised scaling is a linear transformation represented by 
a 2 x 2 matrix that has two real eigenvectors that are not scalar 
multiples of one another. 


3 Diagonalising matrices 


For example, a matrix with two distinct real eigenvalues represents a 
generalised scaling because any eigenvector corresponding to one 
eigenvalue is not a scalar multiple of any eigenvector corresponding to the 
other eigenvalue. Also, dilations are generalised scalings because every 
non-zero vector is an eigenvector of a matrix representing a dilation. 


In contrast, a horizontal shear with non-zero shear factor is not a 
generalised scaling because all eigenvectors of the matrix representing this 


shear are scalar multiples of a 


Activity 17 ‘Investigating generalised scalings g 
Open the Visualising generalised scalings applet. 


Set the first eigenvector to be C) and the second to be G with 


corresponding eigenvalues 1 and 2, respectively. Check that the diagram 
obtained matches Figure 17 (although the colours may differ). 


Change the eigenvectors and eigenvalues to investigate the behaviour of 
generalised scalings. In particular, investigate the effects of positive, 
negative and zero values of the eigenvalues. 


3 Diagonalising matrices 


By now you’ve probably come to see that diagonal matrices have a number 
of properties that make them easy to work with. For example, the 
eigenvalues of a 2 x 2 diagonal matrix are the elements on the leading 
diagonal, and 


(o)  () 


are eigenvectors of the matrix corresponding to these eigenvalues. Also, 
multiplying two diagonal matrices is straightforward, and results in 
another diagonal matrix. 


Here you'll learn about a procedure called diagonalising a matrix, which 
allows you to represent a matrix that has two distinct real eigenvalues in 
terms of a diagonal matrix, and thereby make use of some of the helpful 
properties of diagonal matrices. You'll see that diagonalising a matrix 
gives you a useful method for finding powers of the matrix. Then at the 
end of the section this method will be applied in modelling the populations 
of certain groups of predators and prey. 
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3.1 Diagonalising a matrix 


Suppose that A is a 2 x 2 matrix that has two distinct real eigenvalues. In 
this subsection you’ll learn how to diagonalise matrices of this type, which 
involves finding a 2 x 2 diagonal matrix D and a 2 x 2 invertible matrix P 
such that A = PDP™!. Later in the section you'll see why writing A in 
this way can help you to perform matrix calculations. 


Let’s demonstrate the method for diagonalising a 2 x 2 matrix using the 
matrix 


1 2 
a=( a) 
which you saw in the examples of Section 1 has eigenvalues —1 and 4, and 


3 
that these are indeed eigenvalues and eigenvectors of A without referring to 
the earlier examples by checking the corresponding eigenvector equations: 


Ae) ea 


Next consider the matrix 


1-9 
p=(i 4): 


whose columns are the two eigenvectors of A that we found. Two 
eigenvectors that are scalar multiples of one another must correspond to 
the same eigenvalue, and we can see that neither of the columns of P is a 
scalar multiple of the other. Now, as you saw in Unit 6, this implies that 
the determinant of P is not zero, so P is invertible. In this case, as 


det P = 1 x 3 — 2 x (-1) = 5, 


: 1 2 F 
corresponding eigenvectors e and ( ), respectively. We can check 


the inverse matrix is 


pia te 2) 18 =2 
~dettP\1 1) 5U 1j’ 


We’ll use the inverse matrix P~! shortly. 
If you multiply A on the right by P, then yov’ll notice that the columns of 
the resulting matrix AP are the vectors A (i) and A G: 


apa (3 )(4 3)=( 2) 


Let’s now introduce another matrix, 


ee 
which is a diagonal matrix with the eigenvalues —1 and 4 of A on the 
leading diagonal. The eigenvalue —1 in the first column of D corresponds 
to the eigenvector E of A in the first column of P, and the 
eigenvalue 4 in the second column of D corresponds to the eigenvector 
2 


@ of A in the second column of P. 


If you multiply P on the right by D, then you'll notice that the columns of 


=I 3 


m=-(19(39-(4 


the resulting matrix PD are the vectors — ( i) and 4 Gs 


As you can see, the matrices AP and PD are equal. This is because the 
columns of the first matrix are, from left to right, 


a(i) and a (3) 


and the columns of the second matrix are, from left to right, 


— E and 4 (5) ; 
and the eigenvector equations for A tell us that the first pair of vectors is 
the same as the second pair of vectors. 
If you multiply each side of the equation 
AP = PD 
on the right by P~!, the inverse of P, then you obtain 
APP = PDP~!. 
Since PP~! =I, the 2 x 2 identity matrix, it follows that 
A = PDP. 
You’ve now seen how to write A as a product of three matrices, one of 


them a diagonal matrix. The same procedure can be applied to 
any 2 x 2 matrix that has two distinct real eigenvalues. 


3 Diagonalising matrices 
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Diagonalising a 2 x 2 matrix 


A 2 x 2 matrix A is said to be diagonalisable if there is 
a 2 x 2 diagonal matrix D and a 2 x 2 invertible matrix P such that 


A = PDP. 


The process of finding matrices D and P that satisfy this equation is 
called diagonalising the matrix A. 


If a 2 x 2 matrix A has only one real eigenvalue, and it’s not already a 
diagonal matrix, then it’s not diagonalisable; that is, you cannot find a 
diagonal matrix D and an invertible matrix P such that A = PDP. 


Strategy: 
To diagonalise a 2 x 2 matrix 


Let A be a 2 x 2 matrix with two distinct real eigenvalues. 
1. Find the eigenvalues A; and 2 of A, and corresponding 


eigenvectors i and Cr respectively. 
yı y2 


2. Form the 2 x 2 matrix 


P= 6 
Yl Y2 
3. Find the inverse P~! of P. 
Write down the equation A = PDP™!, where 


_ Pa o 
eG sl 


It’s important to remember that the order in which you place the 
eigenvectors to form the matrix P, and the order in which you place the 
eigenvalues to form the matrix D, must correspond. That is, the first 
column of P must be the eigenvector of A corresponding to the eigenvalue 
in the first column of D (the top-left element), and the second column of P 
must be the eigenvector of A corresponding to the eigenvalue in the second 
column of D (the bottom-right element). If you switch the elements on the 
leading diagonal of D, then you must switch the columns of P (and then 
recalculate P71). 


You should also remember that each eigenvalue of a matrix has many 
corresponding eigenvectors. You can choose any of these eigenvectors to 
place in the columns of P. Just make sure that the two columns of P are 
eigenvectors corresponding to different eigenvalues. 


3 Diagonalising matrices 


| oll 
Example 9 Diagonalising a 2 x 2 matrix A 


Diagonalise the matrix 


oS 


Solution 
@. First find the eigenvalues of A. ® 
We have 
tr A = (—1)+5 = 4, 
det A = (—1) x 5 — 4 x (—2) = 3. 


So the characteristic equation of A is 

A — 4A +3 =0. 
Hence 

(A-1)(A- 3) =9, 
so the eigenvalues of A are 1 and 3. 
@. Find eigenvectors corresponding to these eigenvalues. ® 
The eigenvector equations have the form 

C3" st) Q)-) 

—2 5-A/}\y}) \O}* 


When à = 1, we obtain 


(<2) G) (0). 


This gives 
—2a2 + 4y = 0, 
=2 Aly = O) 


Hence 2x = Ay, so) xr = 2y. ty = 1, then r = 2. Hence o is an 
eigenvector of A corresponding to the eigenvalue 1. 


@. There’s no need to check this now as yov’ll carry out a check at 
the end. & 


When A = 3, the eigenvector equation is 


G 2) G)=(0), 


This gives 
—4z + 4y = 0, 
Hie <= 20) = 0). 
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This pair of equations is equivalent to the single equation 


—x+y=0; thatis, y= zx. 


If c = 1, then y = 1. Hence (i 


i) is an eigenvector of A corresponding 


to the eigenvalue 3. 
@. Form a matrix whose columns are the two eigenvectors. ©& 
Let 
De wl! 
p= 6 7 l 
® Find the inverse of P. ® 


Then 
det RZD me E 


1 1 -1 1 -1 
=i a 
P = detP & 5) e T 


@. Write down the equation A = PDP™!, where D is a diagonal 
matrix with the eigenvalues of A on the leading diagonal. Make sure 
that the eigenvalue of A in the first column of D corresponds to the 
eigenvector of A in the first column of P, and similarly for the second 
columns. © 


We have 
A = PDP}, 


SO 


where 


10 
De ui l 
@. Check that the equation A = PDP! is indeed satisfied. To do 


this, multiply the matrices in the expression PDP, and check that 
you obtain A. ® 


(Check: 


IE DCE Gy 


We could have switched the order of the eigenvalues and eigenvectors in 
the preceding example to give different choices for D and P, namely 


1 2 3 0 
Bei 1 and p=(j a 
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In the activity that follows, and in other activities that require you to 
diagonalise a matrix, you may find that you order the eigenvalues in the 
opposite way to that given in the solution. You may also find that you 
choose different eigenvectors to those given in the solution, affecting the 
matrix P and its inverse. To be sure that the matrices you’ve obtained are 
correct, always check that the equation A = PDP~! is satisfied at the end 
of your working. 


Activity 18 Diagonalising 2 x 2 matrices 


Diagonalise the following matrices. 


@ A=(5 $) wa=(7 3) & 


Taco 
> 
Il 
e 
= N 
l 

aon 
Si 


Jacobi eigenvalue algorithm 


Carl Jacobi was a German mathematician whose research covered 
many mathematical subjects, including calculus, dynamics and 
number theory. One of his numerous contributions to mathematics, 
which came in 1846 while he was working on celestial mechanics, was 
to develop an algorithm to calculate the eigenvalues and eigenvectors 
of a particular type of matrix. 


The method was rediscovered in 1949 by the mathematicians Herman 
Goldstine (1913-2004), Francis Murray (1911-1996) and John von 
Neumann (1903-1957). In his book The Computer from Pascal to von 
Neumann, Goldstine describes how he presented the method at a 
conference in 1951, only to be told that Jacobi had discovered 

it 100 years earlier! 


3.2 Powers of matrices 


One of the uses of diagonalisation is to calculate powers of matrices. To 
see how this is done, let’s begin by investigating powers of diagonal 
matrices, then move on to looking at powers of diagonalisable matrices. 


Suppose that you wish to work out powers of the diagonal matrix 


-(59) 


3 Diagonalising matrices 


Carl Jacobi (1804-51) 
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You will find that 


m= (5 2) (0 = 2): 
m= (5) (0 #)=(0 2) 
m= (5 a) (o a)=(0 2): 


and so on. More generally, we have the following observation. 


Powers of 2 x 2 diagonal matrices 


Let 
a 0 
Dal, a 
Then 
n 
mel Hs tor 7) = 2, Bo ccc 


Activity 19 Finding powers of 2 x 2 diagonal matrices 


Find the following powers of matrices. 


m6 4) oy oo) 


You’ve just seen that there is a straightforward method for calculating 
powers of diagonal matrices; now we’ll use that method to help us find 


powers of diagonalisable matrices. Suppose then that A is 


a 2 x 2 diagonalisable matrix. This means that there is a 2 x 2 diagonal 


matrix D and a 2 x 2 invertible matrix P such that 
A= PDP. 
It follows that 
A? = (PDP!) (PDP t) = PD(P'P)DP™!. 
Since PEP =I, the 2 x 2 identity matrix, we obtain 
A’? = PDIDP ~! = PD’°P '!. 
Similarly, we have 
A? = (PDP!) (PDP H (PDP!) = PDP}, 


and in general the following observation holds. 
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Powers of 2 x 2 diagonalisable matrices 
Let A be a 2 x 2 diagonalisable matrix that satisfies 
A=PDP", 


where D is a 2 x 2 diagonal matrix and P is a 2 x 2 invertible matrix. 
Then 


A” =PD"P! forn=1,2,3,.... 


This gives us a strategy to find powers of diagonalisable matrices. 


Strategy: 
To find powers of 2 x 2 diagonalisable matrices 


Suppose that you want to find A”, where A is a 2 x 2 diagonalisable 

matrix and n is a positive integer. 

1. Diagonalise A; that is, find a diagonal matrix D = ( ) and 
a 2 x 2 invertible matrix P such that A = PDP™!. 


ao o 0 
2. Calculate D = | 


3. Apply the formula A” = PD” P 1. 


Remember that you can only apply this strategy to diagonalisable 
matrices, so that excludes those 2 x 2 matrices that are not diagonal 
matrices and have only one real eigenvalue. 


When you apply the strategy, the final calculation to find PD"P~! may 
involve large numbers, so you might want to check your answer at the end. 
One way to do this is to use the module computer algebra system. Here we 
give an alternative method to check your solution, by hand. This method 
is based on a useful identity about the trace of the product of two 
matrices, which you’ll establish in the next activity. 


Activity 20 Establishing a trace identity 


Given any pair of matrices 


a=(° and Bel? ae 
c d r s 


show that tr (AB) = tr (BA). 


3 Diagonalising matrices 


147 


Unit 11 Eigenvalues and eigenvectors 


Suppose now that P is an invertible 2 x 2 matrix. If you replace A by P7! 
and B by PB in the trace identity tr (AB) = tr (BA), then you obtain 


tr (P~'(PB)) = tr((PB)P~’). 


Since P~'P =I, the identity matrix, the matrix expression on the left of 
this equation simplifies to give B. So we obtain the following useful 
formula. 


Trace formula for 2 x 2 matrices 
If B is a 2 x 2 matrix, and P is an invertible 2 x 2 matrix, then 


iB = i (PBR). 


We can apply this formula to our strategy for working out matrix powers, 
as follows. Suppose that A is a 2 x 2 diagonalisable matrix, so that there 
is a diagonal matrix D and an invertible matrix P such that 


A” =PD"P"! 


for any positive integer n. If we substitute D” for B in the formula from 
the box above, and let the matrix P in the box be the same as the 
matrix P used to diagonalise A, then we obtain 


tr (D”) = tr(PD"P~). 


But PD”P™! is equal to A”, so we conclude that A” and D” have the 
same trace. This gives us a handy way of checking matrix powers, 
summarised below. 


Trace test for checking 2 x 2 matrix powers 


Suppose that the 2 x 2 matrices A, D and P satisfy A = PDP™!, 
where D is a diagonal matrix. Then 


tr (A”) = tr (D”) 


for any positive integer n. 


For example, you saw in Example 9 that if you diagonalise the matrix 


=p 4 
B 


then you obtain the diagonal matrix 
1 0 
D= l l 
Now, you can check that 


aoo —25 52 3 1 0 
< = a Ra a a F 
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Then 
tr(A®) = (-25)4+53=28 and tr(D*)=1427=28, 
so the traces of A? and D? agree, as the trace test predicts. 


Of course, checking that A? and D? have the same trace doesn’t ensure 
that A? = PD?P~!, but nonetheless this trace test is a useful way of 
picking up numerical errors. 


Here’s an example in which the strategy for finding powers 
of 2 x 2 diagonalisable matrices is applied, and the trace test is carried out 
at the end. 


Example 10 Finding a power of a 2 x 2 diagonalisable matrix 


il il 


Diagonalise the matrix A = ( 7 28 


Ib and hence find A®. 


Solution 
@. Diagonalise A. To do this, first find the eigenvalues of A. ® 


Since A is a triangular matrix, the eigenvalues are the elements on 
the leading diagonal, namely —2 and 1. 


@. Find eigenvectors corresponding to these eigenvalues. ©& 


The vector (3) is an eigenvector of A corresponding to the 
eigenvalue 1. 
The eigenvector equations have the form 

Co" 244) =() 

0 —2—A]} \y Ove 

When A = —2, we obtain 

C a) G) =W 

0 O} ly op 

This gives 

aw i Y = 0), 

Ove =- Oy = 0) 
This pair of equations is equivalent to the single equation 


3nr+y=0; thatis, y=—3z. 


If x = 1, then y = —3. Hence E is an eigenvector of A 


corresponding to the eigenvalue —2. 


3 Diagonalising matrices 


@® 
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@. Find a diagonal matrix D and an invertible matrix P such 
that A= PDE # 


Let 
il i 

r=(} 

Then det P = —3, so 
1 -3 -l 
=l = 2 

E ( 0 i) z 
We have 

A = PDP}, 


where 


o=(} 9) 


@. After so many calculations, it’s worth checking that the 
equation A = PDP™! is satisfied, even though you'll carry out 
another check at the end. ® 


(Check: 


(3G Belé a) 


®. Calculate DË. & 
Then 


oa (8 )-( a) 


@®. Apply the formula AÊ = PD°P-!. @ 


So 
= ile 15 (¢ i) 
0 -3/\0 64/3\0 -1 
1 
=3 
1 
3 


@. Check that the traces of A® and DË are equal. ® 
(Check: tr (A) = 1 + 64 = 65 and tr (D®) = 1+ 64 = 65.) 
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Activity 21 Finding powers of 2 x 2 diagonalisable matrices 


—1 0 
2 1 
=2 —6 
2 5 


(a) Diagonalise the matrix A = ( | and hence find A’. 


(b) Diagonalise the matrix A = ( i and hence find A5. 


You may have found that it’s quicker to work out the matrix powers in the 
preceding activity by multiplying the matrices directly rather than using 
the diagonalisation strategy. The benefits in speed for calculating powers 
of matrices using diagonalisation really only become apparent for much 
larger powers. This is because far more computations are needed to work 
out A” than to work out D” when n is a large positive integer, and you 
are more likely to make errors in the former calculation than in the latter. 


Another point to consider here is accuracy. Each time a calculation is 
carried out by a scientific calculator, the accuracy of the answer is limited 
by the number of digits stored in the memory. The more calculations that 
are done, the more errors creep in due to rounding. Calculation of matrix 
powers using diagonalisation involves fewer calculations and hence gives 
more accurate answers. 


3.3 A predator-prey system 


In this subsection you’ll see an application of diagonalisation in modelling 
a wildlife population. Although the model is unrealistically simplistic, and 
the numbers used are made-up, it nevertheless exhibits many of the 
features of more sophisticated models employed by biologists. 


On the Canadian tundra there is a population of Arctic wolves and a 
population of Arctic hares. The wolves rely on the hares to form a major 
part of their diet, so the populations of wolves and hares are related. 


Starting from some fixed time, the populations of wolves and hares are 
recorded after each year passes. Let’s write Wp for the number of wolves 
after exactly n years have passed, and hpn for the numbers of hares after 
exactly n years. In particular, wo and ho are the numbers of wolves and 
hares, respectively, at the very start of the process. We’ll assume that the 
initial number of hares is at least twice the initial number of wolves; that 
is, ho > 2wo. 


An Arctic wolf 


A biologist studying these animals proposes the following system of 
equations to model the populations: 


Wn = 0.94Wy_1 + 0.02hy_1 


24105... 
ha = —0.03wn—1 Ai 4 n ) 
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Each of the coefficients on the right-hand side of these equations indicates 
how much the population of wolves or hares after n — 1 years will affect 
the population after n years. 


For example, the terms 0.94w,,_; and 0.02h,,_1 in the first equation tell us 
that there must be a large number of hares relative to wolves in order for 
the wolf population to grow. Similarly, the terms —0.03w,_1 and 1.01h,_, 
in the second equation tell us that if there are a large number of wolves 
relative to hares, then the hare population will shrink. 


We can represent the system of equations in matrix form as 


wn\  / 0.94 0.02\ [wn E 
(i) > R oD (a) (n =1,2,3,...). 


Let’s write 
a=( oos 10) ad = (72): 
then the matrix equations tell us that 
Xn = Axa (m= l2; 3p): 
Now, as X,_1 = AXņn—2 provided n > 2, we obtain 
Xn = AXn_1 = A(Axn_2) = A?x,,_9. 
Reasoning in a similar way, using X,_2 = Ax,—3 and so on, we find that 
Xn = A’x,_9 


3 
=A Xn-3 
4 
=A Xn—4 


= A”xo. 


To summarise, we have the following formula for calculating the 
populations of wolves and hares after n years. 


Formula for populations of wolves and hares 


The populations of wolves and hares after n years (given by Xn) is 
related to the initial populations of wolves and hares (given by xo) 
according to the formula 


Xn = A”xo. 


Let’s use this formula to predict the populations of wolves and hares 
after 10 years starting from some initial populations that we’ll specify 
later. The first task is to find A!°, and you can do this using the strategy 
of the previous subsection. Here are the results of applying that strategy, 
in brief. You should find that 1 and 0.95 are eigenvalues of A, and that 


1\. ; i , l 
(5) is an eigenvector of A corresponding to the eigenvalue 1, and (7) is 


an eigenvector of A corresponding to the eigenvalue 0.95. Of course, you 
can check these eigenvalues and eigenvectors using the eigenvector 
equations. Next, let 


1 0 is 
D=(j a me P=(5 a 


Since det P = —5, the inverse of P is 


pin If 1 -2\_/-02 04 
Smaps df. 06 —09)* 


Finally, diagonalising A we obtain A = PDP™! and 
A! = PDIP! 


where 
pe- 0 yof D 
~ (o 0.951?) \o0 0.95!) * 


We could now work out the elements of the matrix A!°, but we can 
actually avoid doing so here (remember that it’s x19 that we wish to 
calculate, not A10). As xj9 = Al°xo, it follows that 


x19 = PDP™ txo, 


so to find xj9 you can start by calculating P~!xo, the result of which is a 
vector. You can then multiply this vector on the left by D!° to give 

another vector, and then multiply this vector on the left by P to give xj. 
The next example works through the details for a particular value of xj. 


3 Diagonalising matrices 
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Now 
po! 1 = —0.2 0.4 1 = 1.8 
5 0.6 —0.2 5 —0.4)° 
@. Next find DIOP! : . To avoid rounding errors, leave the 


number 0.951? (the bottom-right element of D) in that form for now, 
and use a calculator to compute numerical values at the end. © 


So 


nome A 0 18) 1.8 
Bap G ~ o 0.95197 (0.4) ~ \—0.4 x 0.9579) ° 
@. Now find PDPP! i & 


Hence 


wN Oo a | 

A E =PD P|. 
{A2 18 
S31) \=04 01952" 


_ {1.8 —-0.8 x 0.957° 
D4 04 0:95.97 


@. Reintroduce the factor 200 that was dropped earlier. ® 


So 
Alx9 = Al? (200 (5) 
_ 10 (1 
- 200” (!) 


1.8 — 0.8 x 0.95" 
ae — 0.4 x nent) 


_ (360 — 160 x 0.9510 
~ \ 1080 — 80 x 0.9519 


@. Work out these values using a calculator. & 
_ (264.20... 
T \ OBO). of 
@. You can’t have 0.2 of a wolf or 0.1 of a hare, so round the final 


answers to the nearest integer. It’s important to round only at the 
end of the calculation. © 


The model predicts that, to the nearest integers, there will 
be 264 wolves and 1032 hares after 10 years. 
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Activity 22 Exploring other initial conditions 


Use the model above to predict the populations of wolves and hares 
after 10 years given that there are initially 600 wolves and 2100 hares. 


Models for predator—prey systems 


The system you’ve seen for modelling populations is rather crude. A 
much more widely used system is provided by the Lotka—Volterra 
equations, which are, 


daz 

dE = x(a = by), 
dy 

a —y(c— dz). 


Here x and y denote the populations of some prey and predator 
species, respectively, at a certain time t (they are functions of t), 
and a, b, c and d are real constants. 


The Lotka—Volterra equations are named after two scientists who 
independently investigated the equations in the early twentieth 
century. The American biophysicist Alfred Lotka first used them in 
the context of chemical reactions in 1910, and the equations were also 
discovered later, in 1926, by the Italian mathematician and physicist 
Vito Volterra, who was carrying out analysis of fish catches in the 
Adriatic Sea. 


Vito Volterra (1860-1940) 


4 Systems of differential equations 


In Unit 8 you learned methods for solving differential equations such as 
dy 
oe 

The general solution of this particular equation is y = Ce’, where C is an 


arbitrary constant. This section is about solving pairs of differential 
equations, such as 


2y. 


dx 

— =I +3 
dt Tt Y, 
dy 

Z = g — 2y. 
d mT 


In this pair of equations, x and y are both functions of t, and the solution 
(to follow later in the section, in Example 13) consists of an expression 
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for x in terms of t and an expression for y in terms of t. You'll see how the 
strategies you’ve learned for finding eigenvalues and eigenvectors can help 
you to solve such differential equations. 


4.1 Matrix differential equations 


A pair of differential equations of the form 


oF arsi 
a Y, 
dy 

q Tt W 


where x and y are functions of a variable t, and a, b, c and d are real 
number parameters, is called a system of differential equations, or 
sometimes just a system. In further mathematics modules, a system of 
differential equations may refer to a more complicated set of equations 
than those considered here (involving, for example, higher derivatives and 
more complex expressions in x and y). 


To save space we will in future write t and y for the derivatives of x and y, 
respectively, with respect to t. With this notation, the system of 
differential equations has the simpler form 

t = ax + by, 

y = cz + dy. 


Using matrix notation we can give another representation of these 
equations. Let 


-0 


Both x and y are functions of t, so x is a vector function of t. We’ll write 
the derivative of x with respect to t as x. It’s given by 


z T 
XS 
(G) 
We can now write the system of differential equations using vectors and 
matrices. 


4 Systems of differential equations 


Matrix form of a system of differential equations 
The system of differential equations 

£=az + by, 

y = cx + dy, 
can be written in matrix notation as 

WC J) 

VY) e gy w 

or more briefly as 

KEA 


where 
A= ( a <= o and x= Ge 
@ @ y y 


Systems of differential equations written in the form x = Ax are called 
matrix differential equations. 


Activity 23 Writing systems of differential equations in matrix 
notation 


Write the following systems of differential equations in matrix notation. 


L=—x+4y t = —y 
(a) Y = Tx — 3y (b) Y= 


A solution of a system of differential equations consists of an 

equation for x in terms of t and an equation for y in terms of t that 
simultaneously satisfy the system. The two equations can be written 
together as a single vector equation. The general solution of the system 
is a solution containing arbitrary constants, such that any solution of the 
system is obtained by choosing particular values for the arbitrary 
constants. Solving a system of differential equations means finding its 
general solution. 


The remainder of this section is concerned with solving systems of 
differential equations. 
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4.2 Decoupled systems 


Consider the system of differential equations, 
t = T, 
y = 4y. 
This is said to be a decoupled system because the first differential equation 


involves the function x but not the function y, and the second differential 
equation involves y but not x. 


The matrix form of this system is 

z\ /—1 0\/z 

v) \ 0 4) W’ 
As you can see, the matrix is a diagonal matrix. In general, you can 
determine whether a system of differential equations is a decoupled system 


by checking whether the matrix in the matrix representation is a diagonal 
matrix. 


Decoupled systems of differential equations 


A decoupled system of differential equations, or more briefly a 
decoupled system, is a system of the form 


ap = (ONE. 
Wy = ay. 


In matrix notation, a decoupled system has the form 
ON fa ON fe 
i 0 a w 


A decoupled system consists of a differential equation in x and a differential 
equation in y, which you can solve separately. Let’s solve the system 
C= 7; 
y = 4y, 
in this way. 
In Unit 8, you learned that differential equations such as t = —x can be 


solved by the method of separating variables. Applying this method 
(with x 4 0) gives 


[ears f-ra 
x 


Now carry out both integrations to obtain 
In |z| = —t + c, 


where c is an arbitrary constant. Hence 


je See: 


4 Systems of differential equations 


It follows that 

a-—ee' (ife>0) and r—-ece! (ifx <0). 
Since e° is a constant, the general solution is 

r= Ce, 


where C is an arbitrary constant. (When C = 0 we obtain the particular 
solution x = 0, which we discounted before applying the method of 
separating variables to avoid division by zero.) 


Similarly, the general solution of the differential equation y = 4y is 
y = De™, 
where D is an arbitrary constant. 


Using similar methods we can find the general solution of any decoupled 
system. 


Solution of a decoupled system of differential equations 


The general solution of the decoupled system of differential equations 


i = ie, 

y = dy, 
is 

r = Ce”, 

y= De®, 


where C and D are arbitrary constants. 


Notice that the coefficients of t in the exponents of e% and e” are equal to 
the diagonal elements of the matrix 


a 0 
0 dj)’ 
Of course, these are also the eigenvalues of the matrix, and this 


observation will be important later, when we consider more general 
systems of differential equations. 
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Example 12 Solving a decoupled system 


Find the general solution of the following system of differential 
equations. 


p= oe, 
Wy = 2. 
Solution 


@. The general solution of a decoupled system of differential 
equations has the form 


r = Ce”, 

y = De*. 
» 
This is a decoupled system of differential equations, so the general 
solution is 

z= Ce’, 

J= De”, 


where C and D are arbitrary constants. 


Activity 24 Solving decoupled systems 


Find the general solutions of the following systems of differential equations. 


(a) 


t = —2r z= 0.22 
j b) 
y = 3y (b) y = —0.5y 


4.3 Coupled systems 


A system of differential equations such as 
&=-—x2+A4y, 
y = —22 + oy, 


is said to be a coupled system because at least one of the two equations 
involves both functions x and y. That is, a coupled system is a system that 
isn’t decoupled. 


4 Systems of differential equations 


The matrix form of this system is 


()-G 


The matrix is not a diagonal matrix, which tells you that the system is a 
coupled system. 


Coupled systems of differential equations 


A coupled system of differential equations, or more briefly a 
coupled system, is a system of the form 


čt = ax + by, 
y = Ce + dy, 


where b and c are not both 0. In matrix notation, a coupled system 
has the form 


()-(¢ YC), 


where ( i) is not a diagonal matrix. 


Let’s discuss how to solve coupled systems. To do this, we recall the 
decoupled system 


LL, 
y = 2y, 


from Example 12. The matrix form of this decoupled system is 
x x 1 0 
(5) =A (5) , where A = ¢ 3) ; 
and we know that 


x, = P is an eigenvector of A corresponding to the eigenvalue A; = 1, 


and 


X= (o) is an eigenvector of A corresponding to the eigenvalue Ag = 2. 


You saw in Example 12 that the general solution to this decoupled system 
is 

r= Cé, 

y= De”, 
where C and D are arbitrary constants. This solution can be written using 
column vectors as 


OOO 
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So you can see that the solution of the decoupled system is of the form 
(5) = Oey + De?*x, 


where x; and x» are eigenvectors of A corresponding to the eigenvalues A, 
and A2 of A, respectively. 


In fact, the general solution of any system of differential equations given by 
a matrix with two distinct real eigenvalues has a similar form. 


Solution of a system given by a matrix with two distinct real 
eigenvalues 


The system 
t = ax + by, 
y=cx+dy, 


d 
and A2, has general solution 


(5) = Cex, + De*?*xo, 


ON PAT s : 
where A = f ) is a matrix with two distinct real eigenvalues \ 1 


y 


where x; and x2 are eigenvectors of A corresponding to the 
eigenvalues A; and Ag of A, respectively, and C and D are arbitrary 
constants. 


Let’s demonstrate why the general solution given in the box above is true, 
using the coupled system 


£t = —x + 4y, 
y = —2x + 5y, 


which was considered at the start of this subsection. The matrix form of 
this system is 


x x —1 4 
@ =A (5) , where A= 2 J . 
In Example 9 you saw that 
©) is an eigenvector of A corresponding to the eigenvalue 1, 


and 
IA : . Í 
,) isan eigenvector of A corresponding to the eigenvalue 3, 


and you saw that you can diagonalise A to give 


A = PDP, 


4 Systems of differential equations 


where 


_ [1 0 _f21 a [1-1 
ba a = i) and P e a 


After diagonalising, the coupled system can be written as 


(i) =P= G) 


Then multiplying both sides of this equation on the left by P~! gives 


p=! a = DP! @ (2) 


Let’s now define 


Then 


D-i G= (234) 


Since x and y are both functions of t, so are u and v. Differentiating both 
sides of this equation with respect to t gives 


(3) = (4355): 


that is 


(2) = (0 36) 


We now have a decoupled system in u and v, rather than a coupled system 
in x and y. You learned earlier how to find the general solution of 
decoupled systems; in this case the general solution is 


U\ C'et — t 1 3t 0 
(3) = (Bas) =o (0) +P (i), 
where C and D are arbitrary constants. 


All that remains is to write the solution for u and v as a solution for x 


and y. Recall that (*) = P7! a which implies that (£) =P J 
v y y v 
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Using the solution for (*) given above, we obtain 


(5) = P(Ce (o) + De” (7) = CtP (o) + De%P (7) 


1 


Since P (o 


) and P (7) are the first and second columns of P, 


respectively, which are the eigenvectors o and (7) of A, we obtain the 


desired solution 


G) -oe (7) roe G) 


Let’s finish by writing this solution in the same form that the system was 
given. To do this, add the vectors on the right-hand side to give 


z\ _ [20e + De’ 
y) (Oe + De% }? 


xr = 2C + De™, 
y = Cé + De*. 


SO 


All this might seem like a lot of hard work, but that’s because we 
explained why the formula for the general solution of this coupled system 
is true; you needn’t go through all these steps each time you wish to solve 
a coupled system. Instead just follow the strategy below. 


Strategy: 
To find the general solution of a system of differential 
equations given by a matrix with two distinct real eigenvalues 


Consider the system of differential equations 

t = az + by, 

ý = cz + dy, 
where a, b, c and d are real numbers such that the matrix J has 
two distinct real eigenvalues. To solve the system, proceed as follows. 


1. Write down the matrix form of the system of differential 
equations. 


2. Find the eigenvalues \1 and A2 of the matrix, and find 
corresponding eigenvectors x; and x2. 


3. Write down the general solution, which is 
—= Ce ox + De™*xs, 
where C and D are arbitrary constants. 


4. State the solution as an equation for x and an equation for y. 


4 Systems of differential equations 


Example 13 Solving a coupled system 


Find the general solution of the following system of differential 
equations 


ie = 2p So oie, 
Q = = = ZY. 
Solution 


@. Write down the matrix form of the system of differential 
equations. ©& 


The matrix form of the system is 


“esal where A= ep 

y yY —] —2 
@. Find eigenvalues and eigenvectors of the matrix. ® 
We have 


tr A = 2 + (—2) = 0, 
det A = 2 x (—2) — 3 x (-1) =-1. 


So the characteristic equation of A is 
X -—1=0. 
Hence A? = 1, so the eigenvalues of A are —1 and 1. 
The eigenvector equations have the form 
E G)=() 
—1 -2-A/)/\y 0} ° 


When A = —1, we obtain 


(4 4) G)=() 


This gives 
3x + 3y = 0, 
a =: 


This pair of equations is equivalent to the single equation 


z+y=0; thatis, y=-—z. 


If x = 1, then y = —1. Hence (i is an eigenvector of A 


corresponding to the eigenvalue —1. 


come AO- 


@ 
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When A = 1, the eigenvector equation is 


(i 3) G)=() 
—1 -3) ly opa 
This gives 
4p oy = O, 
—x¢ — 3y = 0. 
This pair of equations is equivalent to the single equation 


x+3y=0; thatis, x» = —3y. 


If y = 1, then z = —3. Hence z is an eigenvector of A 


1 
corresponding to the eigenvalue 1. 


a(i A)-C) 


@. Write down the general solution. ® 


The general solution is 


ese ob .{-3 
<= Ce (i) + De ( a 


where C and D are arbitrary constants. 


@. Give the solution in the same form that the system was given. In 
this case, write it as two equations. ® 


Hence 
w= Ce IDe 
y= Ge Me: 


Activity 25 Solving coupled systems 


Find the general solutions of the following systems of differential equations. 


t = —2r + 2y t = Tx + 3y 
(a) y = —2xr + 3y (b) y=3r=y 


You’ve now seen how to solve systems of differential equations given by 
matrices with two distinct real eigenvalues. In further mathematics 
modules you can learn how to solve systems given by matrices that don’t 
have two distinct real eigenvalues, as well as how to solve more complex 
systems such as those involving second derivatives. 
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Learning outcomes 


After studying this unit, you should be able to: 


give a geometric interpretation of eigenvalues and eigenvectors 
find eigenvalues and eigenvectors of 2 x 2 matrices 


find eigenvalues and eigenvectors of matrices larger than 2 x 2 using 
the computer algebra system 


use efficient strategies to find eigenvalues and eigenvectors of special 
types of matrices such as triangular matrices and matrices 
representing flattenings, rotations and reflections 


understand the geometric interpretation of special types of matrices 
such as matrices representing shears, flattenings, rotations, reflections 
and generalised scalings 


diagonalise 2 x 2 matrices that have two distinct real eigenvalues 
find powers of matrices by diagonalising 


write systems of differential equations as matrix differential equations, 
and solve those equations given by matrices with two distinct real 
eigenvalues. 
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Solutions to activities 


Solution to Activity 1 
We have 
(e s) G) = (a) =) 
0 3 k 3k k 
Hence 
Ax = 3x. 


Solution to Activity 2 
(a) We have 
(i o) (i) =C) 
1 0/\k]) ~~ kay 
Hence 


Ax = x. 


(b) We have 


G o) (4) = Ca) =a) 


Hence 


Ax = —x. 


Solution to Activity 3 
We have 


—3 -1 1\ | 1\ | 1 1 
4 2) \-4}) ~ \-4}) ~~ \—-4)? 
1\. ; ; 
so i) is an eigenvector of the matrix 


corresponding to the eigenvalue 1. 
Also, 


(a 4) (4) = (2) e) 


1 
corresponding to the eigenvalue —2. 


1\. f ; 
so (_ is an eigenvector of the matrix 


Solution to Activity 4 


(a) The characteristic equation is 


A? — ((—2) + 0)A + (2) x 0 — 1 x (—1) SO. 


That is, 
A? +2 +1=0. 


(b) The characteristic equation is 


X? — ((—3) + (-1))A + (-3) x (-1)-1 x 4=0. 


That is, 
E EEE) 
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Solution to Activity 5 
(a) We have 


(b) 


tr A = (—1) + (—2) = -3, 
det A = (—1) x (-2) -1x0=2. 
So the characteristic equation of A is 

M4 3A+2=0. 
We have 
trA=4+ 
detA=4x5-1xl= 
So the characteristic equation of A is 
NW —8\+1=0. 


Solution to Activity 6 


(a) 


(c) 


We have 
tr A = (-1)+0=-1, 
det A = (—1) x 0 — 2 x 1 = —2. 
So the characteristic equation of A is 
X +A- FM: 
Hence 
(A + 2)(A — 1) = 90, 
so the eigenvalues of A are —2 and 1. 
(Check: (—2) +1 = -1 = tr A.) 
We have 
tr A = (—5) + 1 = —4, 
det A = (—5) x 1 — (-3) x 3 = 4. 
So the characteristic equation of A is 


A? +4A+4=0. 
Hence 
(A +2)? =0, 


so A has a single, repeated eigenvalue —2. 
(Check: (—2) + (—2) = —4 = tr A.) 
We have 
trA=141=2, 
det A=1x1-(-1)x1=2. 
So the characteristic equation of A is 
X — 20 4+2=0. 


Hence 

y- 2) (—2)?-4x1x2 
E 2x1 
_2+y—4 

J 

242i 

=a 

=1+i. 


So the eigenvalues of A are 1 — į and 1 + i. 
(Check: (1-7) +(1+7%) =2=trA.) 
Solution to Activity 7 


The eigenvector equations have the form 


CR 223) G@) = (0): 


When A = —2, we obtain 


G )G)-): 


This gives 
3x + 3y = 0, 
4x + 4y = 0. 


This pair of equations is equivalent to the single 
equation 


x+y=0; thatis, y= —z7. 


1 
eigenvector of the matrix corresponding to the 
eigenvalue —2. 


(Check: G 3 E = (3 =e (i) 
When A = 5, the eigenvector equation is 
Ci -00 
4 —3)}) \y Qy 
This gives 
—4z + 3y = 0, 
4x — 3y = 0. 


This pair of equations is equivalent to the single 
equation 


If x = 1, then y = —1. Hence (i) is an 


3y = 4x. 


If x = 3, then 3y = 4 x 3, so y = 4. Hence e is 


an eigenvector of the matrix corresponding to the 
eigenvalue 5. 


TORORO 


Solutions to activities 


Solution to Activity 8 
(a) We have 


trA=6+2=8, 
det A=6x2-1x5=7. 
So the characteristic equation of A is 
A? -8r+7=0. 
Hence 
(A-1)(A-7) =0, 
so the eigenvalues of A are 1 and 7. 
The eigenvector equations have the form 
C3") Q-@) 
5 2— AJ \y 0J’ 


When A = 1, we obtain 


(N00 


This gives 
ox +y = 0, 
5a +y = 0. 


This pair of equations is equivalent to the single 
equation 


y= =r: 
If x = 1, then y = —5. Hence (_5) is an 


eigenvector of A corresponding to the 
eigenvalue 1. 


cme (3)(8) 


When A = 7, the eigenvector equation is 


Cs -3 G) =): 


This gives 
=p y= 0, 
5x — Sy = 0. 


This pair of equations is equivalent to the single 
equation 


x—-y=0; thatis, y=a. 
If x = 1, then y = 1. Hence C) is an 


eigenvector of A corresponding to the 
eigenvalue 7. 


cine 90-0 
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(b) We have (c) We have 
trA=1+4+3=4, Ax = —4x 
dettA=1x3-0x3=3. for every vector x. So A has only a single, 
So the characteristic equation of A is repeated eigenvalue —4, and any non-zero 
A? —4A+3=0. vector, such as i) is an eigenvector of A 
Hence corresponding to this eigenvalue. 
(A— 1)(A— 3) =9, (d) We have 
so the eigenvalues of A are 1 and 3. trA=2+4+0=2, 
The eigenvector equations have the form det A = 2 x 0—1 x (-1)=1. 
C -=À 3 ) C) = (o) , So the characteristic equation of A is 
oon) AN - —24+1=0. 


When A = 1, we obtain 
Hence 


(0 2) G)=(0)- A-1} =0, 


so A has only a single, repeated eigenvalue 1. 


This gives 
3y = 0, The eigenvector equations have the form 
l ; : : l ' -1 O-A/) \y) \O/}° 
This pair of equations is equivalent to the single a ; 
equation Substituting A = 1 gives 
1 1 x 0 
= 0. = 
(i -1) = (0). 
Choose x = 1. Then (o) is an eigenvector of A This gives 
corresponding to the eigenvalue 1. zr+y=Q0, 
t 3\/1 1 1 —t£—y=0. 
heck: = =1 . 
ee e (o) (o) (5) ) This pair of equations is equivalent to the single 
When \ = 3, the eigenvector equation is equation 
& ) (e) (*) xt+ty=0; thatis, y=-—z. 
E MAY i If x = 1, then y = —1. Hence f is an 
This gives ; 3 
eigenvector of A corresponding to the 
—22%+ 3y =0, : 
eigenvalue 1. 
ee ge | 21\(1\_(1 1 
This pair of equations is equivalent to the single (Check: 1 Gy hep 1 aif -) 
equation 
—2z + 3y = 0; thatis, 3y = 2z. Solution to Activity 9 
We have 
If x = 3, then 3y = 2 x 3, so y = 2. Hence 0 11 1 9 1 
is an eigenvector of A corresponding to the -1 21 Ty] 12) =2 (2) 5 
eigenvalue 3. —1 -1 4 1 2 1 
(Check: € a (3) = (6) =i (3) 0 1 1 1 3 1 
0 6 =1 9 Lt =(3}=3/ 2], 
—1 -l 2 6 2 


170 


0 1 1 2 2 2 
-1 2 1 ipsis 
—1 -1 4 1 il il 
1 1 2 
1], 1 and 1 
1 2 1 


are eigenvectors of the matrix, and the 
corresponding eigenvalues are 2, 3 and 1, 
respectively. 


Solution to Activity 11 


(a) 


Since the matrix is a triangular matrix, the 
eigenvalues are the elements on the leading 
diagonal, namely 1 and 5. 
The vector a is an eigenvector of the matrix 
corresponding to the eigenvalue 1. 
The eigenvector equations have the form 

5-—A 0 x\ /0 

2 1—-A) \y] \O/}° 


When A = 5, we obtain 


G -3 6) = 0) 


This gives 
Oz + Oy = 0, 
2x — 4y = 0. 


This pair of equations is equivalent to the single 
equation 
2x — 4y = 0; 


that is, x = 2y. 


2\. 
1 1S an 
eigenvector of the matrix corresponding to the 


eigenvalue 5. 
aa ( 9) C- (9) -20 


Since the matrix is a triangular matrix, the 
eigenvalues are the elements on the leading 
diagonal, so there is only a single, repeated 
eigenvalue —2. 


If y = 1, then gz = 2. Hence 


1\. : : 
The vector ( i) is an eigenvector of the matrix 


corresponding to the eigenvalue —2. 


(c) 


Solutions to activities 


Since the matrix is a triangular matrix, the 
eigenvalues are the elements on the leading 
diagonal, namely 1 and 3. 
1\. ; i 
The vector (a) is an eigenvector of the matrix 
corresponding to the eigenvalue 3. 
The eigenvector equations have the form 
3-— 1 Ti {0 
0 1-A/ wJ 07’ 


When A = 1, we obtain 


EDG- 


This gives 
2x+y=0, 
Or + Oy = 0. 


This pair of equations is equivalent to the single 
equation 
22+ y = 0; 


that is, y=—2z2. 


If x = 1, then y = —2. Hence (2) is an 


eigenvector of the matrix corresponding to the 
eigenvalue 1. 


om (i D- 


Since the matrix is a diagonal matrix, the 
eigenvalues are the elements on the leading 


diagonal, namely —7 and 7. The vector @ is 


an eigenvector of the matrix corresponding to 


the eigenvalue —7, and G is an eigenvector of 


0 
the matrix corresponding to the eigenvalue 7. 


Solution to Activity 12 

Since the matrix is a triangular matrix, the 
eigenvalues are the elements on the leading diagonal, 
so there is only a single, repeated eigenvalue 1. 


The vector k 


i is an eigenvector of the matrix 


corresponding to the eigenvalue 1. 
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Solution to Activity 13 
(a) We have 


tr A = (—1)+9 = 8, 
det A = (—1) x 9 — (—3) x 3 = 0. 
So the characteristic equation of A is 
X? — 8A =0. 
Hence 
A(A — 8) = 0, 
so the eigenvalues of A are 0 and 8. 


The eigenvector equations have the form 


C's * 3) G) =) 


When A = 0, we obtain 


C3) @)=()- 


This gives 
=g — 3y =0, 
3x + 9y = 0. 


This pair of equations is equivalent to the single 
equation 
x + 3y = 0; 


that is, x = —3y. 


If y = 1, then z = —3. Hence a is an 


eigenvector of A corresponding to the 
eigenvalue 0. 


or (4 =) C) -0-0 


When A = 8, the eigenvector equation is 

(3 DG) =W 

3 1) \y 0} ° 

This gives 

—9x — 3y = 0, 

3x +y = 0. 
This pair of equations is equivalent to the single 
equation 


3xa+y=0; thatis, y= —3z. 


If x = 1, then y = —3. Hence o is an 


eigenvector of A corresponding to the 
eigenvalue 8. 


(Check: 


(3 3) = a) =el 


(Here is an alternative way to find the 
eigenvalues and one of the eigenvectors. 

Since det A = 0, the matrix represents a 
flattening, so the eigenvalues are 0 and tr A, 
and here tr A = 8. Next, since A has 
determinant zero, the first column of A, namely 


is an eigenvector of A corresponding to the 
eigenvalue 8. This eigenvector is not the same 
as the eigenvector 


(3) 


found using the previous method, but both 
vectors are acceptable answers: they are scalar 
multiples of each other using the scale 

factor —1.) 


We have 
trA=14+1=2, 
dettA=1x1-1x1=0. 
So the characteristic equation of A is 
X 2 =0. 
Hence 
A(A — 2) = 0, 
so the eigenvalues of A are 0 and 2. 
The eigenvector equations have the form 
CP )Q-@) 
1 1—A) ly 0} 


When A = 0, we obtain 


(1) G)=()- 


This gives 
r+y=0, 
zt+ty=0. 


This pair of equations is equivalent to the single 
equation 
x+y =O0; 


that is, y= —a. 


1\. 
_;]} isan 
eigenvector of A corresponding to the 
eigenvalue 0. 


(Check: G i) E (o) =0 (4). 
When A = 2, the eigenvector equation is 
Ci) G)-) 
1 -l1/} \y O}° 
This gives 


—-z+y=0, 
xr—y=0. 


If x = 1, then y = —1. Hence 


This pair of equations is equivalent to the single 


equation 


x—-y=0; thatis, y= zx. 


If x = 1, then y = 1. Hence is an 


1 
¢ 
eigenvector of A corresponding to the 
eigenvalue 2. 


a(i 9-6- 


(Alternatively, you can find the eigenvalues and 
one of the eigenvectors using a similar strategy 


to that shown at the end of the solution to 
part (a).) 


Solution to Activity 14 
We have 
trA=0+0=0, 
det A =0 x 0-(-1) x 1=1. 
So the characteristic equation of A is 
AX +1=0. 
Hence 
X? = -1, 
so the eigenvalues of A are —i and i. 
(Check: (—i) +i =0=trA.) 


Solution to Activity 15 


(a) This matrix represents a reflection in a line 
through the origin because it is of the form 


cos 2a sin 2a 
sin2a@ —cos2a)’ 


with a= 7/4. 


eH 


Solutions to activities 


Therefore it has eigenvalues —1 and 1, and 
corresponding eigenvectors 


e = a) 


and 


cos(7/4)\ _ (1/V2 
sin(m/4)) — \1/V/2)”’ 
respectively. Multiplying each of these 


eigenvectors by V2 gives the simpler pair of 
eigenvectors 


(7) = G) 


ome € )()=()--C) 
aa (=) 2) 


This matrix represents a reflection in a line 
through the origin because it is of the form 


cos 2a sin 2@ 

sin2a@ —cos2a}’ 
with a = 7/6. Therefore it has eigenvalues —1 
and 1, and corresponding eigenvectors 


(Tatar) = (3/2) 


and 


eono 


respectively. Multiplying each of these 
eigenvectors by 2 gives the simpler pair of 
eigenvectors 


(a) ma (4) 


cori 
(hs 8) @)-Ca)--@) 
(vain a (i) = (a) =)» 
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Unit 11 Eigenvalues and eigenvectors 


Solution to Activity 16 


(a) Since 


-a)l 


So f maps the point (4,4) to the point (8,8). 


(b) Since 


(D-0 


So f maps the point (—2,—1) to the 
point (—2, —1). 


(c) Since 


So f maps the point (1, —2) to the 
point (—4, —7). 


Solution to Activity 17 


For each eigenvector of the matrix, the 


transformation scales distances in the direction of 
the eigenvector by a factor equal to the magnitude 


of the corresponding eigenvalue. 


If each eigenvalue of the matrix is non-zero, then 


the image of the grid of parallelograms is another 
grid of parallelograms. Each line in the original grid 


is parallel to its image line. 


If the matrix has one positive eigenvalue and one 


negative eigenvalue, then the transformation 


reverses orientation. If the eigenvalues of the matrix 


are both positive or both negative, then the 
transformation preserves orientation. 
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If one eigenvalue of the matrix is zero, and the 
transformation is not the zero transformation, then 
the image grid collapses onto the line through the 
origin in the direction of a vector given by any 
non-zero column of the matrix. In this case the 
transformation is a flattening. 


Solution to Activity 18 


(a) Since A is a triangular matrix, the eigenvalues 


are the elements on the leading diagonal, 
namely —1 and 2. 


The vector (5) is an eigenvector of A 
corresponding to the eigenvalue 2. 


The eigenvector equations have the form 


Po" aha) (= (0). 


When A = —1, we obtain 


EIG- 


This gives 
3x + by = 0, 
Ox + Oy = 0. 


This pair of equations is equivalent to the single 
equation 


3¢+6y=0; thatis, «= —2y. 


1 
eigenvector of A corresponding to the 
eigenvalue —1. 


Let 


=(} 


Then 
det P = 1 x 1 — (—2) x0=1, 


SO 


1 /1 2 1 2 
—1 _ Z= 
as é i) G i 


We have 
A = PDP}, 


where 


If y = 1, then z = —2. Hence (3) is an 


We have 
trA=1+4+2=83, 

dettA=1x2-2x1=0. 

So the characteristic equation of A is 
X —3\=0. 

Hence 
A(A — 3) = 0, 

so the eigenvalues of A are 0 and 3. 


The eigenvector equations have the form 


CD 223) G) =(0): 


When A = 0, we obtain 


(7 2) (7) = (o) 


This gives 

x+2y=0, 

x+2y=0. 
Hence x = —2y. If y= 1, then z = —2. 
Hence ~ is an eigenvector of A 


corresponding to the eigenvalue 0. 


When A = 3, the eigenvector equation is 


Ci =) G)=() 
1 -1/) \y 0j 
This gives 
—2r + 2y = 0, 
g=y =Ù. 
This pair of equations is equivalent to the single 
equation 
that is, 


xz—y=0; Y= 


Iy 
,) is an 
eigenvector of A corresponding to the 
eigenvalue 3. 


If x = 1, then y = 1. Hence 


Solutions to activities 


SO 


We have 
A = PDP}, 


where 


(Alternatively, you can find the eigenvalues and 
one of the eigenvectors of A by using the 
observations about eigenvalues and eigenvectors 
of matrices with determinant zero from 
Subsection 2.2. In this case, as det A = 0, the 
eigenvalues of A are 0 and tr A = 3. 


Furthermore, the vector i) given by the first 


1 
column of A is an eigenvector of A 
corresponding to the eigenvalue 3.) 


We have 
tr A = 2+ (-6) = —4, 
det A = 2 x (—6) — 7 x (—1) = —5. 
So the characteristic equation of A is 
N+ 44-5=0. 
Hence 
(A+ 5)(A — 1) = 90, 
so the eigenvalues of A are —5 and 1. 
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Unit 11 Eigenvalues and eigenvectors 
The eigenvector equations have the form 


Ca a) G)=(0): 


When A = —5, we obtain 


(1-1) () =): 


This gives 

Tx + Ty =0, 

-gy =y =0: 
This pair of equations is equivalent to the single 
equation 

xr+y=0; thatis, y= -—z. 


If x = 1, then y = —1. Hence a is an 
eigenvector of A corresponding to the 
eigenvalue —5. 
When A = 1, the eigenvector equation is 
1 T\) fay) {0 
—1 -7) \y/  \0/’ 

This gives 

z+7y=0, 

—xz —Ty=0. 
This pair of equations is equivalent to the single 
equation 


x+7y=0; thatis, «= —Ty. 


If y = 1, then z = —7. Hence E is an 


eigenvector of A corresponding to the 
eigenvalue 1. 


Let 
1 -7 
p=(2-2). 
Then 
det P = 1 x 1 — (—7) x (-1) = -6, 
so 
Pt =a; i=- Le 
detP \1 1 6\1 1 
We have 
A = PDP}, 
where 


Y 
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Solution to Activity 19 
o (o i) = (e ca) = Co o) 
© (a) (0 2! 

© (o o) =C% o= o 


Solution to Activity 20 


We have 
a b\/p q\ _ fap+br aq+bs 
aB=(° i) (? ae cq+ds 
and 
_(p a\(a b\ _ fpat+qe pb+qd 
Ba=(? J (C Ja a: 
So 
tr (AB) = ap + br + cq + ds 
and 


tr (BA) = pa + qc + rb + sd. 
The expressions on the right-hand sides of these two 
equations are equal, so tr (AB) = tr (BA). 
Solution to Activity 21 
(a) Since A is a triangular matrix, the eigenvalues 
are the elements on the leading diagonal, 
namely —1 and 1. 
The vector (1) is an eigenvector of A 
corresponding to the eigenvalue 1. 


The eigenvector equations have the form 


Ce AG) 


When à = —1, we obtain 


G 2) G=): 


This gives 
Oz + Oy = 0, 
2” + 2y = 0. 


This pair of equations is equivalent to the single 
equation 


2r +2y=0; thatis, y=-—a. 


If x = 1, then y = —1. Hence e is an 


eigenvector of A corresponding to the 
eigenvalue —1. 


Let 


TEN 


Then det P = 1, so 


1 1 0 
Pt = = 
det P G i) 


(1 1): 


We have 
A = PDP}, 
where 
—1 0 
p- (49) 
(Check 


o ((-1)? 0\_/~ı1 0 
= -( o P) iO dy 
That is, D? = D, so 

A? = PD?P™! = PDP™! = A. 


Solutions to activities 


(Yov’ll notice that we chose the matrix P to be 


1 0 0 1 
E ae rather than & 4) 


(with the columns interchanged). This is 
because the left-hand matrix is a lower 
triangular matrix, as is the matrix A. 
Multiplying two lower triangular matrices gives 
another lower triangular matrix, and the inverse 
of a lower triangular matrix is also a lower 
triangular matrix. So choosing P to be a lower 
triangular matrix simplifies the computations in 
the solution.) 
We have 
tr A = (—2) +5 = 3, 

det A = (—2) x 5 — (-6) x 2 = 2. 
So the characteristic equation of A is 

A? —3A+2=0. 
Hence 

(A — 1)(A— 2) = 0, 
so the eigenvalues of A are 1 and 2. 


The eigenvector equations have the form 


Cr” 53) @)=() 
2 5— AJ W 0j- 
When A = 1, we obtain 
CORO 
2 4} \y O}° 
This gives 
—3xr — ôy = 0, 
2r + 4y = 0. 


This pair of equations is equivalent to the single 
equation 


x+2y=0; thatis, z= —2y. 


If y = 1, then z = —2. Hence (1) is an 


eigenvector of A corresponding to the 
eigenvalue 1. 


When A = 2, the eigenvector equation is 


(2 3) @)= 0): 
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Unit 11 Eigenvalues and eigenvectors 


This gives 
—4x — ôy = 0, 
2x + 3y = 0. 


This pair of equations is equivalent to the single 


equation 
2r +3y=0; thatis, 38y = —2z. 
If x = 3, then 3y = —2 x 3, so y = —2. 


Hence e is an eigenvector of A 


corresponding to the eigenvalue 2. 
Let 
—2 3 
P= ( l 2 | 
Then det P = 1, so 
pole 1 —2 -3\__ f2 3 
detP \-1 -2 1 2J” 
We have 
A = PDP}, 


where 


eop 


Then 
15 0 1 0 
— = 
is oes 
So 
AŽ = PDřP™! 


(i 2) 0 CG 2) 
-3 263) 
Ti. 

(Check: tr (A5) 

and tr (D5) = 1 


= (—92) + 125 = 33 
+32 = 33.) 
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Solution to Activity 22 


We have xo = (100) = 300 o 


Now 
pafa [-02 0.4) (2) _ (2.4 
7) \ 06 -0.2)/ \7) \-0.2)' 

So 
‘joey A 0 2.4 
n (7) ae 0.95) @: 


B 2.4 
= -0.2 x 0.9510 / * 
Hence 


10 /2\ _ iop=1 (2 

A (; = PDUP 7 
fi 2 2.4 
NST \—0.2 x 0.951? 


_ (2.4—0.4 x 0.9519 
T (7.2 — 0.2 x 095 7° 
So 


Axo = A! (s00 (;)) 
_ 10 (2 
= 300A G) 


2.4 — 0.4 x 0.9510 
i — 0.2 x es) 


_ (720 — 120 x 0.95'° 
~ \ 2160 — 60 x 0.9519 


_ (648.15... 
= 2124.07...) ` 


The model predicts that, to the nearest integers, 


there will be 648 wolves and 2124 hares 
after 10 years. 


Solution to Activity 23 
» -CI 
o =G 6) 


Solution to Activity 24 
(a) This is a decoupled system of differential 
equations, so the general solution is 
g = Ce“, 
y = De*, 
where C and D are arbitrary constants. 
(b) This is a decoupled system of differential 
equations, so the general solution is 
r= get, 
y= De ~°*t 


where C and D are arbitrary constants. 


Solution to Activity 25 
(a) The matrix form of the system is 
c) =A (;) , where A= e J ; 
We have 
tr A = (—2)+3 = 1, 
det A = (—2) x 3 — 2 x (—2) = -2. 
So the characteristic equation of A is 
X -rA-2=0. 
Hence 
(A+ 1)(A — 2) =0, 
so the eigenvalues of A are —1 and 2. 


The eigenvector equations have the form 


22" 22) ()= (0) 

—2 3-A/\y/ \O/° 
When A = —1, we obtain 

= )G)-@) 
—2 4} \y O}° 

This gives 

—x+2y =0, 

—2r + 4y = 0. 
This pair of equations is equivalent to the single 
equation 

—xz + 2y =0; thatis, x= 2y. 


If y = 1, then z = 2. Hence : is an 


eigenvector of A corresponding to the 
eigenvalue —1. 


(Check: © 5) (7) = (<1) =- (7). 


Solutions to activities 


When A = 2, the eigenvector equation is 


(2 1) G)= (0). 


This gives 
—4¢ + 2y = 0, 
—2r +y =0. 


This pair of equations is equivalent to the single 
equation 


—2r +y =0; thatis, y= 2r. 
If x = 1, then y = 2. Hence >) is an 


eigenvector of A corresponding to the 
eigenvalue 2. 


(Check: C s) @ = a =e Ge 


The general solution is 


_o,-t (2 afl 
x = Ce (7) +06 G: 


where C and D are arbitrary constants. Hence 
x = 20e™ + De”, 
y = Ce™ + 2De”. 


The matrix form of the system is 


(5) =A (C) , where A= é +) ; 

We have 
trA = 7+ (—1) =6, 

det A =7 x (-1)-3x 3 = —16. 
So the characteristic equation of A is 

A? — 6-16 =0. 
Hence 

(A + 2)(A — 8) = 0, 
so the eigenvalues of A are —2 and 8. 


The eigenvector equations have the form 


(a a4) ()=() 

3 —1—- àJ \y 0} ° 
When à = —2, we obtain 

G DEC) 

3 1} \y 0}° 

This gives 

9x + 3y = 0, 

oz +y = 0. 
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Unit 11 Eigenvalues and eigenvectors 


This pair of equations is equivalent to the single 
equation 


3x+y=0; thatis, y=—3z. 
If x = 1, then y = —3. Hence E is an 


eigenvector of A corresponding to the 
eigenvalue —2. 


a(i A) -C-i 


When A = 8, the eigenvector equation is 


(a =) G) = 0) 


This gives 

=z + 3y = 0, 

3x — 9y = 0. 
This pair of equations is equivalent to the single 
equation 


—z + 3y = 0; thatis, x= 3y. 
If y = 1, then z = 3. Hence i) is an 


eigenvector of A corresponding to the 
eigenvalue 8. 


7 8\ 3 24 3 
(Check: C 4) @ = ( s) =8 (7) ) 
The general solution is 


aaf 1 st (3 
x = Ce (3) + De G : 


where C and D are arbitrary constants. Hence 
£ = Ce~** + 3De*, 
y=—38Ce-* + De™, 
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